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Abstract 

Recently, it has become apparent that the thermal stability of topologically ordered systems at finite temperature, 
as discussed in condensed matter physics, can be studied by addressing the feasibility of self-correcting quantum 
memory, as discussed in quantum information science. Here, with this correspondence in mind, we propose a model 
of quantum codes that may cover a large class of physically realizable quantum memory. The model is supported by 
a certain class of gapped spin Hamiltonians, called stabilizer Hamiltonians, with translation symmetries and a small 
number of ground states that does not grow with the system size. We show that the model does not work as self- 
correcting quantum memory due to a certain topological constraint on geometric shapes of its logical operators. This 
quantum coding theoretical result implies that systems covered or approximated by the model cannot have thermally 
stable topological order, meaning that no system can be stable against both thermal fluctuations and local perturbations 
simultaneously in two and three spatial dimensions. 

Keywords: self-correcting quantum memory, thermal stability of topological order, quantum coding theory, 
topological phase, stabilizer formalism, topological quantum field theory 
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1. Introduction 

In recent years, ideas from quantum information science have become increasingly useful in condensed matter 
physics where quantum information theoretical viewpoints give an important insight on studies of many-body entan- 
glement arising in ground states or quasi-particle excitations of correlated spin systems ITT OTI . In particular, it has 
been realized that many interesting physical systems in condensed matter physics may be described in the language of 
quantum codes such as stabilizer codes and its generalization 1 12- 16 1. The emerging closeness between two fields has 
made it possible to study several problems concerning many-body correlated spin systems through quantum coding 
theoretical tools |17|. 

In the present paper, we shall explore such a capacity of quantum coding theory for studying problems in con- 
densed matter physics further, by making use of the correspondence between the following two open problems; 

(a) Feasibility of self-correcting quantum memory. 

(b) Thermal stability of topological order. 

Feasibility of self-correcting quantum memory is an important open problem in quantum information science 
concerning reliable storage of qubits. Thermal stability of topological order at finite temperature is an open problem of 
fundamental importance in condensed matter physics which concerns whether topological order may survive at finite 
temperature or not. While these two problems may look very different from each other, they are fundamentally akin 
to each other, in a sense that by searching for self-correcting quantum memory, one can search for topological ordered 
spin systems which are stable at finite temperature ifTsHSSll . This surprising correspondence enable us to address a 
condensed matter theoretical question, thermal stability of topological order, by analyzing a quantum information 
theoretical question, feasibility of self-correcting quantum memory. 

The main contributions of the present paper toward these two open problems are the followings; 

• We propose a certain model of quantum codes which may cover a large class of physically realizable quantum 
codes, which is called a stabilizer code with translation and scale symmetries. 

• We solve the D-dimensional model exactly by determining its coding properties completely and analyze the 
feasibility of self-correcting quantum memory for D - 1, 2, 3. 

• We establish the connection between self-correcting quantum memory and stable topological order at finite 
temperature, and analyze thermal stability of topological order arising in the model. 
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Figure 1 : The correspondence between the feasibihty of self-correcting quantum memory and the thermal stability of topological order, and the 
plan of the paper. 



The plan of the paper is schematically summarized in Fig. [T] We begin the introduction of the paper by reminding 
why these two problems are important in each field. Then, we give brief summaries of our results in the rest of the 
introduction. 

Feasibility of self-correcting quantum memory: Quantum entanglement decays easily. This underlying diffi- 
culty in quantum information science gave birth to the beautiful art of protecting qubits from decoherence; quantum 
coding theory. The central idea of quantum error- correcting codes is to encode a qubit in many-body entangled 
states, and perform error-corrections so that encoded qubits are not lost. After discoveries of first examples of quan- 
tum codes Ii,12, ,23ri26 | which culminated in stabilizer codes |27|, a large number of quantum codes have been found. 
Now, quantum coding theory constitutes one of the most important building blocks for realizing fault-tolerant quantum 
computation f2^. 

Yet, there still remain important gaps between theoretical constructions of quantum codes and their physical 
realizations as quantum memory devices. First of all, most quantum coding schemes encode qubits dynamically by 
applying a large number of logical gates, and encoding is discussed only in terms of the Hilbert space. However, since 
encoded qubits will be eventually lost in the presence of interactions with the external environment, it is desirable to 
store logical qubits statically in some physical subspaces which are naturally protected from decoherence. 

One plausible approach may be to store logical qubits in the gapped ground space of some quantum many-body 
system. In this light, stabilizer codes constructed with geometrically local generators are promising candidates for 
physical realizations of quantum codes since such local stabilizer codes can be realized as the ground space of gapped 
Hamiltonians by using their local generators as interaction terms. Indeed, from purely theoretical viewpoint, it is 
known that sufficiently accurate and frequent error-corrections can store logical qubits reliably 1 181. However, one 
may still dream of having a quantum memory device which would work without active error-corrections, given the 
difficulties and inefficiencies of performing fast and accurate error-corrections in realitjQ 

Self-correcting quantum memory is an ideal memory device which corrects errors by itself lfT6l [18142111291 . Due to 
the large energy barrier separating degenerate ground states, natural thermal dissipation processes restore the system 
into the original encoded states by correcting errors automatically without any active error-correction. If such a mem- 
ory device could exist, it will be a perfect quantum information storage device which may be used commercially in 
the future. Also, the reliable storage of qubits seems to be the starting point for building scalable quantum computers. 



'To the best of our knowledge, there have been no convincing experimental demonstration of error-corrections. All the demonstrations are 
limited only to phase-type errors, and do not coiTect bit-type errors. Also, it seems very difficult to reach the frequency threshold for reliable 
storage of logical qubits |,18i Finally, it is extremely inefficient to keep performing error-corrections during the time one is storing a qubit if one 
stores a qubit for a long duration such as months or years. 
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There has been significant progress toward construction of self-correcting quantum memory. It has been pointed 
out that the Toric code defined on a four-dimensional system (D = 4) serves as self-correcting quantum memory ifTSl 
l30l . Yet, such a four-dimensional code cannot be embedded in a three-dimensional system. While there have been 
several proposals for three-dimensional self-correcting memory |fT6l |3T1 [32]| . validities of none of these proposals 
have been verified ye|^ In addition, it has been shown that a self-correcting local stabilizer code cannot exist in 
two-dimensional systems EOl |33]| . It seems that two-dimensional spin systems are not promising as resource for 
self-correcting quantum memory. Now, the question we want to address can be summarized as follows: 

• Is it possible to have three-dimensional self-correcting quantum memory? 

Thermal stability of topological order: The feasibility of self-correcting quantum memory is closely related to 
another important open problem in condensed matter physics; the thermal stability of topological order. 

Studies of topologically ordered systems l71[8l [T2l[34H39l have been frontiers of researches in condensed matter 
physics community, as systems with topological order are beyond the description of the Landau's symmetry-breaking 
paradigm which was once considered almost as "theory of everything" for studies of many -body systems. Topologi- 
cally ordered systems are also of practical importance in quantum information science since many-body entanglement 
arising in ground states and quasi-particle excitations of topologically ordered systems is a primary resource for real- 
izing various quantum information processing tasks ll8l fT2l . 

The notion of topological order was originally introduced in order to characterize the stability of ground states of 
many-body quantum systems against local perturbations |35|. Loosely speaking, a system is said to have topological 
order when its ground state properties do not change significantly under any types of small, but finite local perturba- 
tions. This stability of ground states against local perturbations is also valuable for quantum information processing 
since topologically ordered spin systems can be used as good quantum codes with macroscopic code distances BOl . 

However, the situation changes completely when one considers the effect of thermal fluctuations on topologically 
ordered systems. In fact, it is known that topological order in a two-dimensional Toric code is not stable at any finite 
temperature which may be quantitatively seen from the fact that topological order parameters such as topological 
entanglement entropy vanish at any non-zero temperature at the thermodynamic limit |41 1. A similar result is obtained 
in a recent numerical work on topological entanglement entropy in a spin liquid model at finite temperature |42|. It 
seems that topological order in a two-dimensional system is not stable at finite temperature according to general 
studies on the ground state properties of two-dimensional frustration-free Hamiltonians EOl l40l . Now, the question 
concerning the stability of topological order can be summarized as follow^ 

• Is there any system with thermally stable topological order, which is stable against both thermal fluctuations 
and local perturbations? 

Main results: Interestingly, this condensed matter theoretical question on the stability of topological order can be 
addressed through quantum coding theory. In fact, it has been pointed out that when topological order in correlated 
spin systems is stable at finite temperature, such a system can be used as self-correcting quantum memory Il20ll43ll44l . 
This correspondence between self-correcting quantum memory and the thermal stability of topological order may be 
better understood by identifying thermal fluctuations as random errors acting on a quantum memory. Therefore, one 
may take a slight liberty and say that a search for self-correcting quantum memory is equivalent to a search for a 
novel quantum phase with topological order which is stable at finite temperature. 

In this paper, with this correspondence in mind, we analyze coding and physical properties of a certain model of 
local stabilizer codes with physically reasonable constraints |17|. The model, which is called Stabilizer code with 
Translation and Scale symmetries (STS model), is constrained to the following physical conditions. 

• Qubits are defined on a D-dimensional square (hypercubic) lattice with periodic boundary conditions. 

• The Hamiltonian consists only of geometrically local interaction terms with translation symmetries. 

• The number of logical qubits does not grow with the system size (scale symmetries). 



-We shall discuss these proposals in detail in section 3.1 

'The definition of topological order at finite temperature will be clearly stated in section 4 
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We show that the model has a certain dimensional duality on geometric shapes of logical operators, as summarized 
in the following informal theorem: 

Theorem (Dimensional Duality). In a D-dimensional STS model (D < 3), m-dimensional and (D - m)-dimensional 
logical operators always form anti -commuting pairs where m is an integer 

Based on this dimensional duality on logical operators, we give answers to open questions (a) and (b), as summa- 
rized below: 

(a) The three-dimensional STS model does not work as self-correcting quantum memory since the energy barrier is 
finite for the encoding with respect to a two-dimensional logical operator. 

(b) Systems covered or approximated by the model cannot have thermally stable topological order, mearring that 
systems cannot be stable against both thermal fluctuations and local perturbations simultaneously in two and 
three spatial dimensions. 

Code distance: Our results on STS models also provide a partial answer to a long-standing open problem in 
quantum coding theory concerning the upper bound on the code distance of local stabilizer codes. The code distance 
(/ is a measure of the robustness of quantum codes against errors, and one of the ultimate goals in quantum coding 
theory is to find a quantum code with a large code distance for a fixed system size (the total number of qubits). 
While an upper bound on the code distance of stabilizer codes is roughly known B3]| . the upper bound for local 
stabilizer codes is currently not known yet. In other words, despite the fact that the stabilizer formalism is a canonical 
framework in quantum coding theory, we still do not know how robust the best local stabilizer code can be ! 

For a long time, it had been believed that the code distance of local stabilizer codes with qubits is upper 
bounded by 0( ^In): d < 0( Va^) at — > oo since all the examples of local stabilizer codes ever found satisfied this 
upper bound lfT2l[T8l . Later, an example of a local stabilizer code whose code distance scales as 0( y/NlogN) was 
found B6]| . While the code distance of this local stabilizer code exceeds the previously believed upper bound 0{ y/N) 
logarithmically, an example of a local stabilizer code whose code distance exceeds 0{ y/N) polynomially has not been 
found yet. Therefore, the code distance of local stabilizer codes seemed to be upper bounded by 0{N^^^) at A^ — > oo 
where e is an arbitrary small positive number, although no analytical result was known on the upper bound. 

It was recently proven that the code distance of local stabilizer codes is upper bounded by (9(L^ ') where L is a 
linear length of the system and D is the spatial dimension (A^ ~ 0(LP)) [20J. While this work does not rule out the 
possibilities for the existence of a local stabilizer code whose code distance exceeds 0(N^-) polynomially, this bound 
was proven to be tight only for D - 1,2. Thus, whether the tight upper bound is 

d < Oim^') 01 d < 0(L°-^) at N ^ oo (1) 

for D > 2 seems to be one of the most important open questions concerning coding properties of local stabilizer codes. 
Our analysis on STS models provides the following tight upper bound on the code distance of STS models: 

• A three-dimensional STS model has a code distance which is tightly upper bounded by 0{L) where L is the 
linear length of the system. Thus, in a three-dimensional system, the upper bound on code distances turns out 
to be more strict than not only O(L^), but also 0{ y/N). 

Organization: The paper is organized as follows. In section |2] we give a brief review of stabilizer codes and 
introduce STS models. In section [3] we discuss the feasibility of self-correcting quantum memory. In section]?] we 
discuss the thermal stability of topological order at finite temperature. In section |5] we present a certain universal 
property of logical operators shared among all the STS models. In section [6] we discuss possible future problems. 
The definition of topological order at finite temperature is discussed further in [Appendix A A possible relevance 



operators is presented in Appendix C and Appendix D 



to topological quantum field theory is discussed in Appendix B The proof of the dimensional duality of logical 



The main part of the paper is self-consistent and accessible to readers without previous knowledge on quantum 
coding theory and topological order However, the proof part is rather technical, and relies heavily on theoretical tools 
developed previously in LI? 47J by the author. In particular, we owe a lot of arguments to lil7J which introduced and 
solved the two-dimensional STS model originally. 
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2. Stabilizer code witli pliysical constraints 



While local stabilizer codes are physically realizable as quantum memory devices in principle, realistic physical 
systems are often constrained to not only the locality of interaction terms, but also various physical symmetries. 
In this section, we give the definition of the STS model which are local stabilizer codes with translation and scale 
symmetries 1 17|. 

In section 2.1 we give a brief review of stabilizer codes. In section 2.2 we describe the definition of the STS 
model. 



2.1. Stabilizer code 

Here, we give a brief review of stabilizer codes which are quantum codes possessing Hamiltonians to support 
logical qubits in the ground space with a finite energy gap |27|. Some notations which will be used throughout this 
paper are also fixed here. Note that we shall use the notations {) for a set and () for a group. 

Stabilizer formalism: The main idea of stabilizer codes is to encode k logical qubits into physical qubits 
(A^ > k) by using a subspace Vs spanned by states \4i) that are invariant under the action of the stabilizer group S: 

Vs = { |«A> e {C'T" : U\iff) = \^), e S]. (2) 

Here, the stabilizer group S is an arbitrary Abelian subgroup of the Pauli group 

S d r ^ [iI,XuZu...,XN,ZN) (3) 

such that -/ i S. The elements in S are called stabilizers. The logical subspace Vs can be realized as the ground 
space of the following Hamiltonian (see Fig.|2]) 

H^-Y,Sj, S^{SuS2,---) (4) 

j 

since the energy eigenvalue is minimized for states satisfying 5 - \iff} for all j. There are k logical qubits encoded 
in Vs where k = N - G(S). Here, G(S) represents the number of independent generators in S. The ground space is 
separated from excited states by a finite energy gap since eigenstates are simultaneously diagonalized with respect to 
eigenvalues +1 of S j. 



energy 



I gap 



3 - logical qubil 



Figure 2: Energy spectrum of a stabilizer Hamiltonian. Logical qubits are encoded in the degenerate ground space. 

Logical operators: In analyzing properties of logical qubits stored in the ground space, operators called logical 
operators play central roles. Logical operators are Pauli operators which commute with the Hamiltonian, but not 
inside the stabilizer group S. Logical operators can be found inside the centralizer group: 

C ^ {[U e r : [U,Sj] = 0, for all ; } ) (5) 
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which is a group of PauU operators commuting with all the stabilizers. Then, a set of logical operators is 



[U e C : ^ I, U i S' 



(6) 



Logical operators may transform encoded qubits since they act non-trivially inside the ground space Vs- 

Equivalence relation: One may introduce an equivalence relation between logical operators by seeing how they 
act inside the ground space. Two logical operators ( and (' are said to be equivalent if and only if { and {' act in the 
same way inside the ground space: 



o a' e s. 



V|<A> e Vs 



(7) 
(8) 



Therefore, logical operators remain equivalent under multiplications of stabilizers. 

Canonical form: It is often convenient to represent a set of 2k independent logical operators in the following 
canonical form 1147 1 : 



A, 

ri. 
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(9) 



Here, {p and are independent logical operators whose commutation relations are {£p, rp) - 0, [(p, r^] = for p + q, 
[tp, tq\ - and [rp, r^] - 0. Thus, only the operators in the same column anti-commute with each other. Note that 
choices of logical operators are not unique. 

Code distance: The code distance is a measure of the robustness of a quantum code, which is quantified by the 
minimal weight of logical operators: 



d = min(w(t/)) where U e h. 



(10) 



Here, w(U) denotes the number of non-trivial Pauli operators constituting U. The code distance corresponds to a 
minimal number of single Pauli errors necessary to destroy an encoded qubit. Roughly speaking, a quantum code 
with a large code distance can securely protect logical qubits. 

Bi-partition: It is often convenient to split the entire system of qubits into two complementary subsets of qubits 
in studying coding properties of stabilizer codes Il48l |49l . Let us recall a useful formula to study geometric shapes 
of logical operators in stabilizer codes through a bi-partition. For a stabilizer code in a bi-partition, the following 
theorem is known to hold |47| (Fig.|3]l. 

Theorem 1 (Bi-partition). For a stabilizer code with k logical qubits, let the number of independent logical operators 
supported by a subset of qubits R be gR. Then, for an arbitrary bi-partition into two complementary subsets of qubits 
R and R, the numbers of logical operators supported by R and R obey the following constraint: 



gR+gR = 2^. 



(11) 




9r + 9r = 2A: 

Figure 3: A bi-partition of a stabilizer code. Each dot represents a qubit. 



This bi-partition theorem is useful for analyzing geometric sizes and geometric shapes of logical operators. For 
example, if we find a region R where there is no logical operator: gn - 0, we immediately know that all the logical 
operators can be supported inside R since ga = 2k ||20J . Thus, one can restrict geometric regions of qubits where 
logical operators are supported. For an extension of the theorem to subsystem codes, see ||50l |5T1 . 



2.2. Stabilizer code with Translation and Scale symmetries 

Here, we describe the definition of STS models, which are local stabilizer codes with translation and scale sym- 
metries. 

(1) Locality of interaction: Physically realistic systems must have only geometrically local interaction terms. To 
introduce the notion of locality to stabilizer codes, we consider a system of qubits defined on a D-dimensional square 
lattice (hypercubic lattice) which consists of N - Li x ■ ■ ■ x Ld qubits where L,„ is the total number of qubits in the m 
direction for m- 1 , ■ ■ ■ ,D. Therefore, qubits are distributed in the physical space with a metric. 

Here, the entire system is separated into a collection of hypercubes which consists of v = vi x - ■ -xvo qubits without 
overlaps by assuming that «,„ = L^/Vm are integer values (see Fig.|4|i. We consider a block of v = vi x ■ ■ ■ x vo qubits 
as the single unit block which constitutes the entire system. In particular, we consider these unit blocks as single 
composite particles with a larger Hilbert space (C^)®'' (Fig. |4|. Thus, the entire system is viewed as a hypercubic 
lattice of «i x ■ ■ ■ x «£> composite particles. 

Now, we assume that interaction terms of the stabilizer Hamiltonian are defined locally: 



(12) 



where S j are supported inside some regions with 2 x ■ ■ ■ x 2 composite particles (Fig.|4ji. (Otherwise, we coarse-grain 
the system). In this paper, instead of qubits, we consider composite particles as the smallest building blocks of the 
system. 



o o o 
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Figure 4: An illustration of the STS model. A two-dimensional example is shown where a unit block of 3 X 2 qubits is considered as a composite 
particle with a larger Hilbert space. Interaction terms S j are defined locally inside a region of 2 X 2 composite particles. The Hamiltonian is 
invariant under unit translations of composite particles. 



(2) Translation symmetries: Physically realistic systems often have not only local interactions, but also some 
physical symmetries. Here, we assume that the stabilizer Hamiltonian possesses translation symmetries: 



T„(H) ^ H (m= I, - - ,D) 
where T„, represent unit translations of composite particles in the m direction (Fig.Hl) 



(13) 
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For simplicity of discussion and in order to accommodate translation symmetries, we set the periodic boundary 
conditions. Then, the entire system may be viewed as a D-dimensional torus: = S' x ■ ■ ■ x where is a circle. 
Thus, the entire system has a topologically non-trivial geometric shape a priori. 

(3) Scale symmetries: In this paper, we are interested in coding properties at the limit where the system size goes 
to the infinity (in other words, at the thermodynamic limit). So far, we have considered the cases where the system 
size n = (ni, - ■ ■ ,no) is fixed. Here, we consider changes of the number of composite particles n,„ while keeping 
interaction terms S j the same. 

It is commonly believed that there is a tradeoff between the number of logical qubits k and the code distance 1* 401 
ISTI where the code distance d decreases as the number of logical qubits k increases for a fixed A^. Then, it may be 
legitimate to limit our considerations to the cases where the number of logical qubits k remains small when the system 
size increases. 

We assume that stabilizer codes have scale symmetries by requiring that the number of logical qubits kpt is inde- 
pendent of the system size n: 

kfi = k, Mn. (14) 

Here, we emphasize that, in a system with scale symmetries, the number of logical qubits k remains constant under 
not only global scale transformations: ft cn where c is some positive integer, but also arbitrary changes of n„,. 

One might think that scale symmetries are too strong as physical constraints. However, through appropriate 
coarse-graining, a fairly large class of local stabilizer codes with translation symmetries can be considered as the STS 
model. For example, let us consider the cases where the number of logical qubits kj^ is small: 

kg < ko, Vit (15) 

where kjj does not grow with the system size n. Then, there always exists some finite coarse-graining such that 
a coarse-grained system has scale symmetries, as proven in |17|. As a result, by analyzing coding properties of 
stabilizer codes in the presence of scale symmetries, one can easily deduce coding properties of stabilizer codes 
which satisfy Eq. ([TSj. Therefore, solutions of STS models are sufficient to discuss coding properties of translation 
symmetric stabilizer codes with a small number of logical qubits. Note that STS models cover the Toric code and its 
generalizations to D-dimensional systems. 

Translation equivalence of logical operators: There is a certain property of logical operators which emerges 
naturally as a result of translation and scale symmetries. For STS models, the following theorem holds ifTTl (Fig.|5]l. 

Theorem 2 (Translation equivalence). For each and every logical operator I in an STS model, a unit translation of 
{ with respect to composite particles in any direction is always equivalent to the original logical operator I: 

T„{{) ~ {, "il e Lg (m = I,-- - ,D) (16) 

where Lg is a set of all the logical operators for an STS model defined with the system size ft. 

Translation equivalence 



-I 



Figure 5: The translation equivalence of logical operators. Each square represents a composite particle. Shaded regions represent translated logical 
operators which are equivalent to each other 

Here, we only give an intuition on why this theorem holds. Let us consider the case where the system size n is 
large. Then, since the number of logical qubits k does not depend on the system size, k is relatively small compared 
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with the system size ft. Now, due to the translation symmetries of the system Hamiltonian, translations of a given 
logical operator £ are also logical operators. However, there are only 2k independent logical operators. Then, there 
must be a finite integer such that { ~ T"n'{() for all the logical operators {. (Otherwise, there would be so many 
independent logical operators). It turns out that a„, - 1 for any { and m. While we have used only the condition that 
the number of logical operators k is small, due to scale symmetries {k is constant), one can prove the above theorem 
by showing a,,, = 1 for any {, m and it. 



3. Feasibility of self-correcting quantum memory 

In this section, we present one of our main results in this paper, concerning coding the feasibility of self-correcting 
quantum memory. In section 3.1 we begin by giving brief discussion on the relation between the qubit storage time 
and the energy barrier. In section 3.2 we describe geometric shapes of logical operators in STS models, and discuss 
whether a three-dimensional STS model works as self-correcting quantum memory or not. 



3.1. Self-correction, energy barrier and storage time 

In subsection, we review how self-correcting quantum memory works by establishing the connection between the 
energy barrier and the qubit storage time. 

Self-correcting classical memory: For simplicity of discussion, let us start by analyzing an example of self- 
correcting classical memory. Consider two-dimensional Ising model: 

H = -^Z,jZ,+ij-^Z,jZ,j+i (17) 
'../ 

which consists of Lx L qubits with periodic boundary conditions. The model works as a classical code since one can 
encode a classical bit in the ground space by labeling |0 ■ ■ ■ 0) as and |1 ■ ■ ■ 1) as 1. 

Now, let us see why this model works as self-correcting classical memory. Assume that the system is originally 
|0 ■ ■ ■ 0). Then, in order for errors to change a ground state |0 ■ ■ ■ 0) into another ground state |1 ■ ■ ■ 1), errors must flip 
all the spins from |0) to However, during these spin flips, the excitation energy becomes at least 0(L) because 
there is a domain wall separating the regions with |0)s and |l)s (Fig.|6]l. In other words, ground states |0 - ■ -0) and 
|1 • • • 1) are separated by a large energy barrier. Then, before errors accumulate, natural thermal dissipation processes 
restore the system into the original encoded state|^ Therefore, the system corrects errors by itself. 

One may estimate the bit storage time of two-dimensional Ising model by using the so-called Arrhenius law: 

T ~ exp(A£'/r) (18) 

up to some polynomial corrections where r is the storage time, AE is an energy barrier and T is the temperature. 
According to the Arrhenius law, the bit storage time can be estimated as r ~ EXP(L) since the energy barrier is 
AE ~ 0(L). While it should be noted that this is an empirical law, it is commonly believed that the law correctly 
estimates the bit or qubit storage time of spin systems (for systems with well defined energy barrier) below a critical 
temperature. Indeed, the law is rigorously verified in various models of classical and quantum memory |fT9l 12111221 
l52l l53ll . However, it should be noted that verification of the Arrhenius law for classical and quantum memory is 
usually a very difficult task, involving an evaluation of time evolution of the system in a presence of interactions with 
an external environment. 

Next, let us consider an example of a classical code which is not self-correcting; one-dimensional Ising model: 

H = -YjZjZj,,. (19) 
j 



Precisely speaking, the system does not return to the original state, but return to a state which is sufficiently close to the original state with a 
probability approaching to unity at the thermodynamic limit. Therefore, one can reliably read out the encoded bit from such a state. 
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One may see that this model does not work as self-correcting classical memory since one can change from |0 ■ ■ ■ 0) to 
1 1 ■ ■ ■ 1) by costing only a finite energy. In other words, one can create a kink by costing only a finite energy: 

10000- ■■0) ^ 10111 ■■■!> ^ 10011 ■■■1) ^ > 11111 ■■■1) (20) 

which may propagate the lattice freely without costing any extra energy. Then, according to the Arrhenius law, the bit 
storage time is t ~ (9(1) which is independent of the system size since the energy barrier is AE ~ (9(1). 
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Figure 6: How self-correcting classical memory works in two-dimensional Ising model. 



Energy barrier and logical operators: One can associate the self-correcting property of two-dimensional Ising 
model with geometric shapes of logical operators by viewing the model as a stabilizer code. Note that D-dimensional 
Ising model satisfies the definition of STS models since interaction terms are local and translation symmetric, and 
there is a single logical qubit {k = 1 and d - 1) regardless of the system size. Two-dimensional Ising model has the 
following pair of zero-dimensional and two-dimensional logical operators: 

e = Zi,i, r = Yl^ij. (21) 
'J 

Then, a classical bit is encoded in eigenstates of a zero-dimensional logical operator £. In order to change the encoded 
bit, one needs to apply a two-dimensional logical operator r since |1 ■ ■ ■ 1) - r|0 • • • 0). Then, an intermediate state 
during the change from |0 ■ ■ ■ 0) to 1 1 ■ ■ ■ 1) may be represented as r*|0 ■ ■ ■ 0) where r* is some "subpart" of the original 
two-dimensional logical operator r (see Fig|7]^a)). Since interaction terms anti-commute with Pauli operators at the 
boundary of r*, the excitation energy associated with r*\Q ■ ■ ■ 0) is proportional to the perimeter of r*. Thus, during 
the change from |0 ■ ■ ■ 0) to |1 ■ ■ ■ 1), the excitation energy must become 0{L) since the perimeters of subparts of a 
two-dimensional logical operator r are always one-dimensional. 

One can also understand why one-dimensional Ising model does not work as self-correcting classical mem- 
ory through geometric shapes of logical operators. One-dimensional Ising model has the following pair of zero- 
dimensional and one-dimensional logical operators: { - Zi and r = 11; ^i- Since a subpart r* of a one-dimensional 
logical operator r is always zero-dimensional, the energy barrier is AE ~ (9(1) (see Fig|7|b)). 

Self-correcting quantum memory: Next, let us discuss how self-correcting quantum memory works. It is 
known that many of good local stabilizer codes with macroscopic code distances do not have self-correcting proper- 
ties since they have string -like logical operators which lead to (9(1) energy barrier As an example, let us consider 
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(b) 



I X : x I X , X X x^ X , x I X ; x"1 



Figure 7: (a) A subpart of a two-dimensional logical operator. The excitation energy is 0(L). (b) A subpart of a one-dimensional logical operator. 
The excitation energy is 0(1). 



two-dimensional Toric code: 



H 



Toric 



(22) 



where qubits live on edges of Lx L square lattice. Interaction terms A, and Bp are shown in Fig. [8] Two-dimensional 
Toric code can be viewed as an STS model since the number of logical qubits is constant: k - 2 regardless of the 
system size. 
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Figure 8: (a) A subpart of a two-dimensional logical operator. The excitation energy is 0(L). (b) A subpart of a one-dimensional logical operator. 
The excitation energy is 0(1). 



One can understand why two-dimensional Toric code does not work as self-correcting quantum memory through 
geometric shapes of logical operators. Since the Toric code has pairs of one-dimensional logical operators as shown 
in Fig. [8] the energy barrier is AE ~ 0(1). As a result, the qubit storage time is t ~ (9(1) according to the Arrhenius 
law. A similar discussion holds for three-dimensional Toric code too. 

As two-dimensional and three-dimensional Toric codes have 0(1) qubit storage time, they are not reliable quantum 
memories. Then, what kinds of quantum codes may serve as reliable quantum memory devices? It is worth noting 
that two-dimensional Toric code works as a reliable quantum memory device if one perform sufficiently frequent and 
accurate error-corrections. In particular, it has been shown [ 18 1 that the qubit storage time can be exponentially long: 
T ~ EXP(L) in the presence of active error-corrections; however, it seems very difficult to reach such accuracy and 
frequency which are necessary for reliably storing qubits. 

In the same paper 1 18 1, it has been also pointed out that four-dimensional Toric code may have exponentially long 
storage time r ~ EXP(L) below the critical temperature since the model has a large energy barrier which scales as 
0(L). This remarkable insight has been further investigated in ||30l . and later, rigorously verified in Il22l . One can 
associate the self-correcting property of four-dimensional Toric code with geometric shapes of logical operators, as 
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the model has pairs of two-dimensional logical operators which lead to 0{L) energy barrier. Yet, since one cannot 
embed four-dimensional Toric code in a three-dimensional space, one hopes to have three-dimensional self-correcting 
quantum memory whose qubit storage time grows as the system size increases. 

Previous proposals: So far, there have been several interesting proposals of three-dimensional self-correcting 
quantum memory. First, a certain model of three-dimensional subsystem codes (16], called three-dimensional quan- 
tum compass model, was proposed as a candidate of self-correcting quantum memory. Subsystem codes may be 
considered as generalizations of stabilizer codes which are supported by Hamiltonians whose interaction terms are 
Pauli operators, but may anti-commute with each other While the model opened a new possibility of quantum mem- 
ory devices supported by frustrated Hamiltonians and initiated studies of Hamiltonian realizations of subsystem codes, 
its validity as self-correcting quantum memory has not been verified since it is difficult to solve the frustrated Hamil- 
tonian, and its properties are hard to determine. Also, it seems difficult for such frustrated systems to have the thermal 
stability since the model supports gapless excitations. Furthermore, its validity has been denied in several literature 
including lfT9l . 

Later, an interesting model of the mixture of two-dimensional Toric code and three-dimensional Bosonic gas, 
called the Toric -Boson model, was proposed 1,31 J where Bosonic gas induces confining potential between anyonic 
excitations. While the model opened new capacities of quantum codes constructed in the so called mixed-dimensional 
configurations, which are of particular interest in ultracold atom physics community, the model itself has two serious 
drawbacks as a candidate of self-correcting quantum memory. First, its storage time is only polynomial in L: t ~ 
POLY(L) since an effective energy barrier is: AE ~ LOG(L). However, it is known that it takes at least 0(d) gate 
operations, which cannot be implemented simultaneously, to read out the encoded qubit, where d is the code distance 
of the system. Then, since the code distance d scales polynomially in L in the model, polynomially long storage time 
is not sufficient for a model to work as efficient quantum memory device ET\ . Second, the model needs a very strong 
coupling between the Toric code and Bosonic gas whose strength scales polynomially in L. Therefore, the model has 
a problem in ffie scalability. 

Recently, an interesting proposal of three-dimensional spin glass models ||32l has been made whose relaxation 
dynamics is very slow due to the existence of a large number of energy local minima. The model seems to have a 
polynomially long storage time: r ~ POLY(L) with logarithmically large energy barrier: AE ~ LOG(L), though it 
needs verifications (so, it could be longer or shorter). The system undergoes a phase transition at T = 0, which mat 
imply a potential thermal instability of the system properties. Finally, it seems difficult o find an efficient decoding 
scheme for the model. Here, it should be noted that the model does not have scale symmetries since the number 
of logical qubits is highly sensitive to the system size, and there is no finite upper bound on it. See section |6] for 
discussion on stabilizer codes without scale symmetries too. 

At this moment, validities of none of these proposals of three-dimensional self-correcting quantum memory has 
been verified yet, and the feasibility of self-correcting quantum memory remains open. Below, we summarize the 
qubit storage time and the energy barrier in these proposals of self-correcting quantum memory. 



Code 


Energy Barrier 


Storage Time 


Comment 


2-dim Toric Code 


0(1) 


0(1) 




3 -dim Toric Code 


0(1) 


0(1) 




2-dim Toric Code 


0(1) 


EXP(L) 


with error-correction 


4-dim Toric Code 


0(L) 


EXP(L) 


self-correcting 


3-dim subsystem code 


? 


7 


gapless excitation? 


Toric-Boson model 


LOG(L) 


POLY(L) 


POLY(L) coupling 


Spin glass model 


LOG(L) 


POLY(L) 


r = phase transition 



3.2. Dimensional duality of logical operators 

In this subsection, we describe geometric shapes of logical operators in STS models, and discuss the feasibility of 
self-correcting quantum memory. 

One-dimension: We begin by reviewing geometric shapes of logical operators in one-dimensional STS models. 
In a one-dimensional chain of ni composite particles, we denote jth composite particle as Pj. Then, the following 
theorem holds [17i . 
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Theorem 3 (Dimensional Duality). For a one-dimensional STS model, there exists a canonical set of logical opera- 
tors: 



{\, ■ ■ ■ , 4 
n, ■ ■ ■ , rk 



(23) 



whose pair of anti-commuting operators (j and rj has one of the following property (j— 1, ■ ■ ■ ,k). 



• £j is a zero-dimensional logical operator defined inside Pi, while rj is a one -dimensional logical operator 
defined in a periodic way: Ti (rj) — rj. 

It is worth presenting geometric shapes of logical operators graphically (Fig. |9]). The code distance of one- 
dimensional models is upper bounded by v which is the inner dimension of a composite particle. Since there always 
exist zero-dimensional logical operators, such systems cannot have topological order. 
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Figure 9: Logical operators in a one-dimensional system. 



Two-dimension: Next, we analyze geometric shapes of logical operators for two-dimensional STS models. Let us 
first introduce some regions of composite particles in order to define geometric shapes of logical operators concisely. 



A square region of xi x X2 composite particles is denoted as P(xi , X2) (Fig. 10 a)): 



P(xi,X2) = [Pr,,r, : 1 < ri < XI, 1 < r2 < X2 } (24) 

where 1 < xi < ni and 1 < X2 < na- Note that a composite particle at the position (ri, r2) is denoted as P,,,^,- Then, 
the following theorem holds II17I . 



Theorem 4 (Dimensional DuaUty). For a two-dimensional STS model, there exists a canonical set of logical opera- 
tors: 



{\, ■ ■ ■ , 4 

r\, ■ ■ ■ , rt 



(25) 



whose pair of anti-commuting operators (j and rj has one of the following two properties (j — 1, ■ ■ ■ ,k). 



£j is a zero-dimensional logical operator defined inside P(l,2v), while rj is a two-dimensional logical operator 
defined in a periodic way: Ti(rj) = rj and 7^2 (o) ~ 

tj is a one -dimensional logical operator defined inside P(l,n2) in a periodic way: 7^2 (4) — 4'. ^f^'l^ ^ 
one-dimensional logical operator defined inside P(«i, I) in a periodic way: Ti(rj) — rj. 



It is worth presenting geometric shapes of logical operators graphically (Fig. 10 b)). There is a dimensional duality 
on geometric shapes of logical operators as follows: 



dim, 1 dim 
2 dim, 1 dim 

14 



(26) 



Logical operators are periodic in the directions in which they are stretched. The code distance is upper bounded by 
0{L) which is consistent with the upper bound OilP^"^) obtained in EOl . The existence of one-dimensional logical 
operators gives rise to topological order. 



(a). 



X2 




(b) 0-dim 



l-dim 



2v 




2-dim 



1-dim 




Figure 10: Dimensional duality in a two-dimensional system, (a) A region of x\ X X2 composite particles is denoted as P{x\,X2). (b) Geometric 
shapes of logical operators in a two-dimensional STS model. 



Three-dimension: Finally, let us proceed to coding properties of a three-dimensional STS model. A region with 
^1 X ^2 X XT, composite particles is denoted as P{xi,X2, X3): 

P{xuX2,Xi) = {P,i,/-2,f3 : 1 JCm, = 1,2,3 } (27) 

where Pri,r,,rj. represents a composite particle at (ri,r2,r3). Then, for logical operators in a three-dimensional STS 
model, the following theorem holds. 

Theorem 5 (Dimensional Duality). For a three-dimensional STS model, there exists a canonical set of logical oper- 
ators: 



■ ■ ■ , 4 
n, ■■■ , n 



(28) 



whose pair of anti-commuting operators (j and rj has one of the following four properties. 

• {j is a zero-dimensional logical operator defined inside P(l,2v, (2v)^), while rj is a three-dimensional logical 
operator defined in a periodic way: T\{rj) — rj, Txifj) - rj and Tj,(rj) - rj. 

• £j is a one-dimensional logical operator defined inside P(ni,2v, I) in a periodic way: T\{{j) — {j, while rj is a 
two-dimensional logical operator defined inside P(\,n2, n^) in a periodic way: T2(rj) — rj and T^(rj) — rj. 

• £j is a one-dimensional logical operator defined inside P(\,n2,2v) in a periodic way: T2((j) — (j, while rj is a 
two-dimensional logical operator defined inside P{ni, 1 , n^) in a periodic way: T\ (rj) — rj and T^irj) — rj. 

• {j is a one-dimensional logical operator defined inside P(2v, 1, ns) in a periodic way: T^{{j) — {j, while rj is a 
two-dimensional logical operator defined inside P(«i, ^2, 1) in a periodic way: T\{rj) — rj and T2(rj) — rj. 

This is the main technical result of the present paper. We present the proof of the theorem in appendices |Appendix| 
[C]and Appendix D It is worth representing geometric shapes of logical operators graphically (Fig. 111. Note that 
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logical operators are periodic in the directions in which they are stretched. There is a dimensional duality on geometric 
shapes of logical operators as follows: 



Odim, 
3 dim, 



1 dim 

2 dim 



(29) 
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1-dim 



1-dim 



l-dim 




Figure 1 1 : Dimensional duality in a three-dimensional system. 

As a result of this dimensional duality, one may find the upper bound on the code distance. When n \ 
the code distance of a three-dimensional STS model is upper bounded as follows: 



d < 2vL 



0{L). 



n2-n-i-L, 
(30) 



Note that this bound is tight for the three-dimensional Toric code. The system may have topological order when it is 
free from zero-dimensional logical operators. 

Higher-dimensions: Though our primary interests are in coding properties of three-dimensional STS models, it is 
worth extending the analysis to higher dimensions. For a D-dimensional STS model {D > 4), we make the following 
conjectures: 

• In a D-dimensional system, m-dimensional and {D - m)-dimensional logical operators form anti-commuting 
pairs where m is an integer 

• The code distance is tightly upper bounded by O(L^) when D is even and by O(L^) when D is odd. 

Note that generalizations of the Toric code to D-dimensional systems have the above dimensional duality for 



arbitrary integer m. See Appendix B.2 for constructions of D-dimensional Toric code. Such a dimensional duality 
arising in geometric shapes of logical operators may be viewed as a manifestation of Poincare duality on coding 



properties of stabilizer codes supported on a D-torus, as very briefly discussed in Appendix B 



However, we report that straightforward generalizations of analysis tools presented in this paper did not work for 
STS models with D > 4, and the conjecture above needs to be verified. 

Feasibility for a three-dimensional STS model: Finally, let us show that three-dimensional STS models do 
not work as self-correcting quantum memory. First, in order for the system to work as a quantum code, there must 
be an anti-commuting pair of one-dimensional and two-dimensional logical operators. We denote them as i and r 
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respectively. Let us assume that the system is initially in the eigenstate of ^ = 1 denoted as - 1)). Then, 
- -1)) - r\t//({! = 1)) and r is a two-dimensional logical operator. Since the boundary of a subpart of r is one- 
dimensional, the excitation energies associated with intermediate states are 0{L). Thus, the encoding with respect to 
( is self-correcting. 

Next, let us assume that the system is initially in the eigenstate of r = 1 denoted as \tff(r = 1)). Then, \t//(r = -1)) - 
(\4i{r - I)) and £ is a two-dimensional logical operator. Since the boundary of a subpart of { is zero-dimensional, 
the excitation energies associated with intermediate states are Oi l). Thus, the encoding with respect to r is not self- 
correcting, and the qubit storage time is r ~ (9(1), independent of the system size. Therefore, while such a system 
is a good quantum code with the code distance 0{L), it works only as self-correcting classical memory, with the bit 
storage time t ~ EXP(L). However, in order for the system to work as self-correcting quantum memory, there must 
be a pair of anti-commuting two-dimensional logical operators. 

Summary and discussion: We summarize coding properties of STS models based on dimensions of pairs of 
logical operators: 



Spatial dim Logical operators Code distance Bit storage time Qubit storage time Self-correction 



1 dim 


dim + 1 dim 


0(1) 


OH) 


0(1) 




2 dim 


dim + 2 dim 


Oil) 


EXP(L) 


0(1) 


classical 


2 dim 


1 dim + 1 dim 


Oil) 


OH) 


0(1) 




3 dim 


dim + 3 dim 


Oil) 


EXP(L) 


0(1) 


classical 


3 dim 


1 dim -t- 2 dim 


Oil) 


EXP(L) 


0(1) 


classical 


4 dim 


2 dim -t- 2 dim 


OiL^) 


EXP(L) 


EXP(L) 


quantum 



where, for D = 4, we presented coding properties of four-dimensional Toric code. 



4. Thermal stability of topological order 

In this section, we present another main result of the present paper, concerning the thermal stability of topological 
order at finite temperature by establishing the connection between the feasibility of self-correcting quantum memory 
and the stability of topological order 



In section 4. 1 we begin by establishing the connection between self-correcting classical memory and the thermal 



stability of ferromagnetic order at finite temperature. In section |4.2| we establish the connection between quantum 



codes and topological order at zero temperature. In section 4.3 we establish the connection between self-correcting 
quantum memory and topological order at finite temperature, and analyze the thermal stability of topological order 
arising in STS models. Note that discussion in this section is rather heuristic, and more rigorous treatment follows 



in Appendix A 



4.1. Classical equivalence 

In this subsection, we establish the connection between self-correcting classical memory and thermal stability of 
ferromagnetic order at finite temperature: 

Classical code <-> Ferromagnetic order aiT -Q 

Self-correcting classical memory Ferromagnetic order at T > 

Two-dimensional ferromagnet: We have seen that two-dimensional Ising model works as self-correcting clas- 
sical memory since the energy barrier separating two degenerate ground states is OiL). This self-correcting property 
of two-dimensional Ising model is closely related to the thermal stability of ferromagnetic order as seen from the 
following thermodynamic argument. 

Consider two-dimensional Ising model: 

Hie) = - 2 - Yj ^'■•^■^y+i " ^ Z ^'''i ^^1^ 

U i,j 
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with an initial bias (symmetry -breaking field); -eYjij^u e > 0. As a result of an initial bias, the system is not 
degenerate anymore, and the ground state is |0 ■ ■ ■ 0). The thermal stability of ferromagnetic order can be analyzed 
through the expectation value of the total magnetization m: 



at finite temperature, which can be computed as follows: 

('«) 



,. 1 51ogZ 
= hm - - 



de 



(32) 



(33) 



where the partition function is: Zip, e) - Tre^^^^'^K Here, we evaluate the expectation value of m at the limit where 
e — > after we take the limit where goes to infinity. (Otherwise, there is no use of introducing an initial bias). Then, 
we have: 



1 > 



'■J ' e^O 

* 'J 



> 



= 



(T = 0) 

{Tc >T > 0) 

(T > T,) 



(34) 
(35) 
(36) 



where Tc is some finite transition temperature, as plotted in Fig. 12 a) 



One may notice that the total magnetization has some non-zero value as long as the temperature is below the 
transition temperature Tc- This indicates that the system is stable against thermal fluctuations for T < Tc- In particular, 
the system properties for Tc > T > are close to the ground state properties at T = 0. However, the system properties 
change drastically at T - Tc, and for T > Tc, the total magnetization vanishes. This indicates that the system 
undergoes a phase transition atT = Tc, and is unstable against thermal fluctuations for T > Tc- 



(a) 2D Ising 
m A 




(b) 1D Ising 
m 



To 
Thermally stable 



Tc = 



1 



♦►T 



Thermally unstable 



Figure 12: Total magnetizations, (a) Two-dimensional Ising model, (b) One-dimensional Ising model. 



Now, one can establish the connection between self-correcting properties and the thermal stability. In a coding 
theoretical language, the computation of the magnetization at finite temperature may be interpreted as follows. We 
first apply an initial bias and encode a bit to a ground state |0 ■ ■ ■ ). Then, we let the system interact with the 
external environment at finite temperature, and wait until the system reaches the equilibrium while slowly turning off 
the initial bias. Finally, we decode the initially encoded bit by measuring the total magnetization. Therefore, self- 
correcting property of two-dimensional Ising model with exponentially long bit storage time is a direct consequence 
of the thermal stability of ferromagnetic order at finite temperature. 

In general, the total magnetization we have used above is called an order parameter in condensed matter physics 
since it can be used to distinguish different phases of matters. On the other hand, quantum coding theoretically, one 
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may view the total magnetization as a symmetric summation of a zero-dimensional logical operator Z,j- with which a 
classical bit is encoded. Later, we shall see that expectation values of logical operators can be used as order parameters 
to characterize the stability of topological order. 

One-dimensional ferromagnet: Next, let us see that one-dimensional Ising model, which does not work as 
self-correcting classical memory, is not stable against thermal fluctuations. Consider one-dimensional Ising model: 



with an initial bias, and compute the expectation values of the total magnetization. Then, we have 



V 2 

j 

V 2 



1 (r = 0) 



= (r > 0) 



(37) 

(38) 
(39) 



as plotted in Fig. 12 b). One may notice that the total magnetization becomes zero at any finite temperature. This 
indicates that the system is not stable against thermal fluctuations at any finite temperature, and the system undergoes 
a phase transition at T = as seen from the sudden change of the magnetization at T = 0. The system properties 
for r = and for T > are significantly diff'erent, and the ground state properties are not stable against thermal 
fluctuations. This is the underlying reason why one-dimensional Ising model has 0(1) bit storage time, and does not 
work as self-correcting classical memory. 

Summary of the equivalence: With observations above, one may notice that a large energy barrier, which is 
essential to self-correcting properties, is the key to the thermal stability of ferromagnetic order. In Fig. 13 we give a 



summary of the equivalence concerning classical memories. However, it is not clear whether a large energy barrier is 



sufficient for the thermal stability of the system or not. We address this issue briefly in Appendix A. 3 
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' Ferromagnetic 



V (Perturbation) 



Figure 13: A summary of the classical equivalence. 



Here, it is worth noting that ferromagnetic systems, such as one-dimensional and two-dimensional Ising model 
are not stable against perturbations, as described in Fig. [13] If one chooses V as the initial bias, one can easily break 
the ground state degeneracy, and thus, the ground state property is not stable against perturbations. This indicates that 



ferromagnetic systems are not topologically ordered (see section 4.2 too). 



4.2. Existence of topological order and quantum codes 

Now, we return to the main discussion in this section, concerning the stability of topological order at finite tem- 
perature. In this subsection and the next, we establish the connection between self -correcting quantum memory and 
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the thermal stability of topological order: 



Quantum code 
Self-correcting quantum memory 



Topological order at T = 
Topological order at T > 



In this subsection, we describe the definition of topological order in spin systems on a lattice, and argue that quantum 
codes have topological order at zero temperature. Note that our discussion closely follows pioneering works B4ll54l . 

Phenomenological definition: Let us first begin by describing the phenomenological definition of topological 
order which is commonly used in physics community. A system is said to have topological order when its ground 
state properties do not change significantly under any types of small perturbations 1351 . 

Definition 1 (Stability against perturbations). Consider a degenerate spin system defined on some closed geometric 
manifold governed by a Hamiltonian H. The system is said to be topologically ordered at T = if and only if the 
system satisfies the following conditions: 



There exists some finite positive number 5 such that, for any perturbations V: 



H' = H+V, 



V 



Vi 



(40) 



where Vj are locally defined and \Vj\ < 6, the ground state degeneracy is not broken at the thermodynamic limit 
(Fig.O. 

The perturbed ground space G can be approximated by the original ground space Go through some local unitary 
transformation: 



uGoW ~ a 



which can be represented as 

U - I exp(-iht)dt 
Jo 

where h is some hermitian operator which is a summation of local terms with finite amplitudes. 



(41) 



(42) 



Energy 



Energy 



|sap 



Original System ; H 



Perturbed System : H +V 



Figure 14: The stability of the ground state properties against local perturbations. 

It is worth noting that a classical ferromagnet (Ising model) is not topologically ordered while it works as self- 
correcting classical memory. First of all, the ground state degeneracy of a classical ferromagnet is broken if a local 
magnetic field term Zj is added as a perturbation: 



V 



(43) 
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While the original Hamiltonian has |0 ■ ■ ■ 0) and |1 ■ ■ ■ 1) as ground states, the perturbed Hamiltonian has |0 ■ ■ ■ 0) as a 
single ground state for positive e. Second, the ground state property significantly changes as a result of perturbations. 
While a classical ferromagnet has -^(|0 ■ ■ ■ 0) + 1 1 ■ ■ ■ 1)) as a ground state, the perturbed ground state is |0 ■ ■ ■ 0), and 

there is no local unitary transformation which transforms :^(|0 ■ ■ ■ 0) + |1 ■ ■ ■ 1)) into |0 ■ ■ ■ 0) 1 1 1. As a straightforward 
extension of this discussion, one may notice that STS models with zero-dimensional logical operators are not topo- 
logically ordered, and topological order may exist only in two or higher-dimensional systems due to the dimensional 
duality of logical operators. 

Coding theoretical definition: While the definition above captures the stability of ground state properties in 
topologically ordered systems, one may wonder what kinds of systems are actually topologically ordered. Also, one 
may wonder if there is a universal feature commonly shared among all the topologically ordered systems. 

Somewhat surprisingly, a quantum coding theoretical viewpoint gives a hint to answer these questions by pro- 
viding an alternative definition of topological order which may capture universal properties of topologically ordered 
systems UJ. In fact, a system of spins on a lattice may be said to be topologically ordered when the ground space 
is separated from excited states by a finite energy gap at the thermodynamic hmit and the ground space realizes a 
quantum code whose code distance d is comparable to the system size f5T|. One may rewrite this characterization of 
topological order through coding properties more explicitly in the following way|^ 

Definition 2 (Quantum code). Consider a degenerate spin system defined on some closed geometric manifold. The 
system is said to be topologically ordered when it satisfies the following condition: 

• Let the degenerate ground state space be G, which is protected by a finite energy gap. Consider reduced density 
matrices of ground states for a region R which is contractable to a single point: 



So, is a topologically trivial zero-dimensional region (Fig. 15 i. Then, all the reduced density matrices of 
degenerate ground states are the same at the thermodynamic limit 

pRim = pRiW)) for all |iA>, e G at N ^ ^. (45) 



Let us see why the condition above ensures that the system serves as a quantum code with a macroscopic code 
distance. When all the reduced density matrices are the same for all the topologically trivial regions, any errors acting 
inside such regions cannot change the encoded logical qubits. In order to change logical qubits, some global operator 
acting on a region with a topologically non-trivial geometric shape is necessary. Such a global operator has a weight 
equal to the code distance. Therefore, the ground space of such a system works as a good quantum code with a 
macroscopic code distance]^ 

Equivalence of two definitions: Now, let us see that definition|2]implies definition[T] Under a small perturbation 
V which is a summation of local terms, we perform a perturbative analysis. Since the ground states are separated by 
excited states by a finite energy gap, only the coupling between degenerate ground states is dominant. Since there is no 
local operator to couple two ground states, one needs to consider (9(L)th order perturbation to have a coupling between 
degenerate ground states. Since such a coupling will be exponentially suppressed by a factor exp(-L/Lo) where Lo 
is some finite length scale, and the energy splitting is exponentially suppressed and goes to zero as L goes to infinity. 
Thus, the ground state degeneracy is protected at the thermodynamic limit. When the ground state degeneracy is 
protected and the ground space is separated from excited states by a finite energy gap, one can approximate the 
ground space from the original unperturbed ground space according to the theory of adiabatic continuation 1541 l55l . 
A complete mathematical proof was presented in lt54J . However, it is not known if definition[T]implies definition|2]or 
not. 



'Strictly speaking, tliis definition is not complete. See 154] for a complete definition. 

*The definition used here is slightly different from the one used in [44j where the existence of the energy gap is not required. The difference 
comes apparent whether a two-dimensional Bacon's subsystem code 1161 is classified as topologically ordered system or not. The ground state 
space of this subsystem code has the same reduced density matrix for any zero-dimensional regions. However, it is not likely to have a finite energy 
gap since it is highly frustrated, and thus, its coding properties are likely to be unstable against perturbations. 
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degenerate ground states reduced density matrices 




Figure 15: A coding theoretical definition of topological order. 



Here, we briefly mention the usage of the expression "topological order" since it is sometimes used in a fuzzy 
and elusive way. In the present paper, by topological order, we mean the ground state properties of topologically 
ordered systems. By its definition, topological order is stable against any types of small local perturbations, meaning 
that the ground state degeneracy is protected and perturbed ground states can be approximated through local unitary 
transformations. Therefore, by the ground state properties, we mean any global features of ground states which are 
free from length scales and are not affected by local unitary transformations. We will return to discussion on the 
meaning of "global features" later. 

Topological order in STS models: Let us conclude this subsection on the definition of topological order by 
seeing that two-dimensional and three-dimensional STS models can have topological order according to definition |2] 
Consider the case when two-dimensional STS model has k pairs of one-dimensional logical operators. Let us split the 
entire system into two complementary regions P{n i - 1 , «2 - 1 ) and P(n i - 1 , «2 - 1 ) where P{n i-l,n2-l)isa region 
with «! - 1 x«2- 1 composite particles. Note that - 1, n2 - 1) is topologically trivial zero-dimensional region since 
one can shrink it into a single point continuously. Since there are 2k logical operators defined inside P(ni - 1, n2 - 1)^ 
we have g f(„^_i n^_i) - 2fc. However, this means that gp(n,-\,n2-i) - from theorem|Tj and there is no logical operator 
defined inside P{ni - 1, n2 - !)■ Then, all the reduced density matrices are the same, and such a two-dimensional STS 
model is topologically ordered since it satisfies definition [2] and thus, definition [T] A similar discussion holds for a 
three-dimensional STS model with pairs of one-dimensional and two-dimensional logical operators. 

4.3. Quantum equivalence 

In this subsection, we establish the connection between self-correcting quantum memory and topological order at 
finite temperature, and analyze the thermal stability of topological order arising in STS models. 

Two-dimensional Toric code: Two-dimensional Toric code is topologically ordered since it has a macroscopic 
code distance: d ~ (9(1). Yet, it does not work as self-correcting quantum memory since the energy barrier is 
AE ~ 0(1). This is closely related to the thermal instability of topological order as we shall see below. Let us perform 



22 



a thermodynamic analysis on two-dimensional Toric code in a way similar to Ising model by adding an initial bias: 



Htie) = Hjoric-eYj^iit) (46) 

L 

HAe) = HToric - eY^T[{r) (47) 



where ( and r are one-dimensional logical operators extending in the 1 and 2 directions respectively. Here, we define 
the following "normalized logical operators": 



L 

1 

mi 



= ^Yj^iit) (48) 
1 ^ 



1=1 

by taking symmetric summations of logical operators, as we did for Ising model. Then, we have 

<m^),^o = <"V>.^o = 1 (r = 0) (50) 
{mt),^o = <m,),^o = (r > 0) (51) 



as plotted in Fig. 16 a) where {m[) is evaluated for Hc{e), and (w,) is evaluated for H,-{e). This indicates that the 
system is not stable against thermal fluctuations at any finite temperature, meaning that topological order arising in 
two-dimensional Toric code is thermally unstable. This implies that one cannot read out initially encoded qubit by 
measuring nif and m,-. 

Three-dimensional Toric code: Next, let us consider thermodynamic properties of three-dimensional Toric code. 
While three-dimensional Toric code does not work as self-correcting quantum memory, it works as self-correcting 
classical memory. This is because it has pairs of one-dimensional and two-dimensional logical operators, and as 
a result, the bit storage time is t ~ EXP(L) while the qubit storage time is t ~ 0{l). These coding properties 
are closely related to thermodynamic properties of three-dimensional Toric code, as seen from expectation values of 
logical operators: 

1 ^ 



where ^ is a two-dimensional logical operator extending in the 1 and 2 directions, while r is a one-dimensional logical 
operator extending in the 3 direction. Their expectation values behave as follows: 

<m,),^o = 1 (T ^ 0) (54) 
<m,),^o = (r > 0) (55) 



and 



<m,),^o = 1 (T ^ 0) (56) 
1 > <m,),^o > (r, > r > 0) (57) 
<m,),^o = (T > T,) (58) 
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where Tc is some finite transition temperature, as plotted in Fig.[T6|b). This imphes that three-dimensional Toric code 
undergoes phase transitions both at T = and T - Tc, and the ground state properties are not completely stable against 
thermal fluctuations at any finite temperature. Yet, the ground state properties partially survive at finite temperature 
as a direct consequence of being self-correcting classical memory. 



(a) 2D Toric 

mi, Mr 



4 



(c) 4D Toric 



Tc = 



To 




(b) 3D Toric 



Tc = 



Too 



iTc 



Figure 16: Expectation values of logical operators, (a) Two-dimensional Toric code, (b) Three-dimensional Toric code, (c) Four-dimensional Toric 
code. 

Four-dimensional Toric code: Finally, let us see that four-dimensional Toric code, which works as self-correcting 
quantum memory, is stable against thermal fluctuations. Here, we define 



1 ^ 



(59) 



(60) 



where ^ is a two-dimensional logical operator extending in the 1 and 2 directions, and r is a two-dimensional logical 
operator extending in the 3 and 4 directions. Then, we have 



{T = 0) 

(r, > r > 0) 

(7^ > T,) 



(61) 
(62) 
(63) 



where Tc is some finite transition temperature, as plotted in Fig. 16 c). This implies that the ground state properties 
are stable against thermal fluctuations and topological order arising in four-dimensional Toric code is stable at finite 
temperature. 

Summary of the equivalence: With these observations, one may notice that large energy barrier inside the ground 
space, which is essential to self-correcting properties, is the key to the thermal stability of topological order. In Fig. 17 
we give a summary of the equivalence concerning quantum memory. 

With this connection between the feasibility of self-correcting quantum memory and the thermal stability of topo- 
logical order, we conclude that topological order arising in STS models is not stable against thermal fluctuations. In 
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Figure 17: The quantum equivalence. 
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other words, for D < 3, there is no system which is stable against both thermal fluctuations and local perturbations 
simultaneously. While discussion here is rather heuristic, we give a more rigorous treatment on the definition of the 
thermal stability of topological order in Appendix A We summarize physical properties of STS models based on 
dimensions of pairs of logical operators: 

Spatial dim Logical operators Local perturbations Thermal fluctuations Memory property 



1 dim 

2 dim 

2 dim 

3 dim 

3 dim 

4 dim 



dim + 1 dim 

dim + 2 dim 

1 dim + 1 dim 

dim + 3 dim 

1 dim + 2 dim 

2 dim + 2 dim 



stable 

stable 
stable 



Classical code 
stable Classical self-correction 

Quantum code 
stable Classical self-correction 

Quantum code 
stable Quantum self-correction 



where, for Z) = 4, we presented coding properties of four-dimensional Toric code. 



5. Emergence of topology in logical operators 

So far, we have addressed specific questions concerning coding and physical properties of gapped spin systems on 
a lattice. While these two questions are of particular importance in quantum information science and condensed matter 
physics, the ultimate goal is to find universal properties of arbitrary gapped spin systems on a lattice and develop a 
unified theoretical framework to discuss their coding and physical properties. In particular, it would be very beautiful 
if arbitrary gapped spin systems can be analyzed universally through some unified theoretical tool which is yet to be 
discovered. 

In this section, as a necessary first step toward this goal, we make an attempt to find uiuversal properties of 
logical operators which are commonly shared among all the STS models. Logical operators in STS models have 
a certain interesting topological property which is a direct consequence of physical constraints we have imposed on 
stabilizer Hamiltonians. Here, we show that one can always deform geometric shapes of logical operators continuously 
while keeping them equivalent by applying appropriate stabilizers as long as one does not change geometric shapes 
in a topologically non-trivial way. We prove this continuous deformability of logical operators in STS models for 
D - 1, 2, 3 with various examples of continuous deformations of logical operators. Some of these results were also 
presented in our previous papers lfT7ll47ll . especially for D -2. 



The program of finding a universal theory of gapped spin systems is continued in Appendix B where the role of 
topology in analyzing coding and physical properties of stabilizer codes is further examined, while we concentrate on 
demonstrating topological properties of logical operators in this section. 

One-dimension: We begin by illustrating a topological property of logical operators for one-dimensional STS 
models. Recall that P{x) represents a region of x composite particles, and gp^i) -km one-dimensional systems 
regardless of the system size since there always exist k zero-dimensional logical operators. Then, as a result of the 
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bi-partition theorem (theorem[T]i, we have 



(64) 




Figure 18: A deformation of a zero-dimensional logical operator. 



Now, let us interpret the equation above from a geometric viewpoint. Consider some logical operator £ defined 
inside P{1). Then, since gp^i) = gpjj^ - k, there must exist some equivalent logical operator {' ~ { which is defined 

inside P{1). In other words, if there exists a logical operator defined inside P(l), one can shrink its geometric shape 
into f (1) by applying some stabilizer (Fig.fTS): 

PW P{1). (65) 



Here, one may notice that P(l) and P(l) are topologically equivalent regions since one can deform P(l) into P(l) 
continuously. On the other hand, if a logical operator is a one-dimensional logical operator defined all over the lattice 
(i.e. defined inside P{n) where n is the linear length of the system), one may not be able to deform the logical operator 
intoF(l). 

Two-dimensions: Next, we analyze a topological property of logical operators for two-dimensional STS models. 
For the convenience of presentation, we assume that zero-dimensional logical operators defined inside P(l,2v) can 
be actually defined inside P(l, 1) in a two-dimensional STS model. This may be done through some appropriate 
coarse-graining. 

We first list regions which serve as references to classify geometric shapes of logical operators in a two-dimensional 
system. We define topological unit regions as follows (Fig.[T9|a)): 



2(0,0) = P( 1,1), 2(1,0) = P(ni,l), 2(0,l) = P(l,n2), Q(l,l) = P(nun2). (66) 

"1" and "0" represent whether a region extends in the corresponding direction or not. For example, Q(l,Q) and 
2(0, 1) are one-dimensional unit regions which extend in the directions of 1 and 2 respectively. 2(0^0) is a zero- 
dimensional unit region with a single composite particle. 2(1' 1) is a two-dimensional unit region which consists of 



all the composite particles in the system. These topological unit regions are shown graphically in Fig. 19 a). We also 
denote a union of all the m-dimensional topological unit regions as R^: 

7?„ = 2(0,0), /?! =2(1,0) U 2(0,1), R2 = Q(l,l). (67) 

We call R„ m-dimensional concatenated topological unit regions. All the concatenated unit regions are graphically 
shown in Fig.[T9|b). 

In a two-dimensional system, there are five different unions of topological unit regions: R^), 2(1,0), 2(0, 1), 
and 7?2. We call these regions, except R2, reference regions, whose set is denoted as 

Rref^{Ro,Q(UO),Q(OA),Ri}. (68) 
Then, one can introduce equivalence relations between these reference regions and their complements in terms of 



continuous deformations. For example, as shown in Fig. 20 a), Rq can be continuously deformed into Ri by enlarging 



the hole of Rq gradually. Also, as shown in Fig. 20 b), Ri can be continuously deformed into Rq since both Ri and 



^0 are zero-dimensional regions without any winding around the torus. Finally, as shown in Fig. 20 c), 2(1,0) can 



be deformed into 2(1,0) since both regions have a winding in the 1 direction. In summary, we have the following 
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Q(i,i) 



Figure 19: Reference regions, (a) Topological unit regions, (b) Concatenated unit regions. 



equivalence relations between reference regions and their complements: 



Ro^Ru Ri^Ro, G(1,0)- 2(1,0), 2(0, 1) - 2(0, 1). 



(69) 



(a) 



□ 
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Figure 20: The topological deformations of logical operators. 



Now, we discuss how geometric shapes of logical operators can be determined. A useful observation regarding 
geometric shapes of logical operators can be obtained by considering the number of independent logical operators 
defined inside a region R. Let the number of independent logical operators inside R be gR. Here, we consider the case 
where we have two regions R and R' where R is larger than R', meaning that R includes all the composite particles 



inside R' (Fig. 21 1. Then, if gR = gR', R and R' support the same logical operators since all the logical operators 
defined inside R have equivalent representations which are supported inside R'. This means that, for a given logical 
operators £ defined inside R, one can always find another equivalent logical operator £' defined inside R'. In other 



words, one can deform the geometric shape of { into f by applying some appropriate stabilizer (Fig. 21 1. Thus, by 
finding two connected regions R and R' where R is larger than R' and g^ - gRi, one can conclude that logical operators 
defined inside R can be deformed into R'. 

With the above observation on geometric shapes of logical operators and deformations in mind, let us describe 
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Qr — 9R' 



Figure 21: A shrinkage from R to R' wlien gn = gR' 



a topological property of geometric shapes of logical operators. The following lemma summarizes the topological 
property of logical operators in two-dimensional STS models. 



Lemma 1 (Topological shrinkage). In two-dimensional STS models, the following equations hold: 



(1,0) 



gQiOA) - 8q(0J) 



k. 



(70) 



In other words, one can shrink geometric shapes of logical operators by applying some appropriate stabilizers in 
the following ways: 



Ri, 



R^ 



Ro, 



2(1,0)^ 2(1,0), 2(0,1)^ 2(0,1). 



(71) 



Note that these shrinkages preserve topological properties of geometric shapes of logical operators. 



Proof. Let us begin with the proof of - gR^ . Let ko be the numbers of pairs of anti-commuting zero-dimensional 
and two-dimensional logical operators. Let ki be the number of pairs of anti-commuting one-dimensional logical 
operators. Then, we notice that gR, - 2k\ + ko since Ri supports both zero-dimensional and one-dimensional logical 
operators. On the other hand, since g^ - 2k - gR^ from theorem |T| and - ko, we have g-^ - gR,- If we use 
theorem|T]to g^ = gR,, we readily obtain gR^ - g^^ too. Next, let us show gg^jy - gQ(i). Note that = k. Then 



we have ; 



k, and have 



□ 



Three-dimensions: Let us continue our analysis on STS models for higher-dimensional cases (D > 2). Below, 
we analyze the topological property of logical operators in three-dimensional STS models. We assume that zero- 
dimensional logical operators and one-dimensional logical operators in a three-dimensional STS model can be defined 
inside P(l, 1, 1), P(«i, 1, 1), P(l,«2, 1) and P(l, 1,«3). 

We begin by finding reference regions for D-dimensional systems. Reference regions for D > 2 can be defined 
from topological unit regions in a way similar to two-dimensional cases. Let d be an arbitrary binary D component 
vector d - (d\, - ■ ■ , do) with dm = 0, 1. Then, topological unit regions are: 



Q(ct) = P(x), where 



(72) 



For example, Q(l, 1,0) = P{n\,n2, 1) for D - 2>. We denote the weight of the binary vector d as w{d) = Y!^t=\ dm, 
which represents the dimension of Q(d). Then, concatenated unit regions are defined as follows: 



Rm= U 

vv(J)=m 



(73) 



A set of reference regions can be obtained by considering all the possible unions of Q{^, which is denoted as Rref- 
One may notice that topological unit regions Qid) are like independent generators for mth homology group for a 
D-torus: Hm(T"). 
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It is worth presenting some examples here. In a three-dimensional system (D - 3), we have the following 
topological unit regions. 



Odim: 2(0,0,0) 

1 dim: e( 1 , 0, 0), Q(0, 1 , 0), 2(0, 0, 1 ) 

2 dim: 2(1,1,0), 2(0,1,1), 2(1,0,1) 

3 dim: 2(1,1,1)- 

Some examples are shown in Fig.|22] Also, concatenated topological unit regions are: 

2(0,0,0) 

Ri =2(i,o,0)u2(0,i,0)u2(0,o,i) 
= 2(1, 1,0) u 2(0, i,i)u 2(1,0,1) 
= 2(1,1,1) 



(74) 



(75) 



which are described in Fig. 23 



One can introduce the equivalence relations in terms of reference regions. Equivalence relations among them are 



29 




Figure 23: Concatenated unit regions in a tliree-dimensional system. 



shown as follows: 

Ro - Ri 

e(i,o,0) - 2(1, 1,0) u 2(1,0,1) 
e(0, 1,0) - 2(1, 1,0) u 2(0, 1,1) 
2(0, 0, 1) - 2(1,0, i)u 2(0, 1,1) 
2(1,0,0) u 2(0, 1,0) - 2(1, 1,0) u 2(0,0,1) 

2(0, 1 , 0) u 2(0, 0, 1) - 2(0, 1 , 1) u 2(1 , 0, 0) (76) 
2(0, 0, 1) u 2(1, 0, 0) - 2(1,0, i)u 2(0, 1,0) 

Ri^Rl 

2(1,1,0)^2(1,1,0) 
2(0, 1,1) -2(0, 1,1) 
2(1,0,1) -2(1,0,1). 

Then, we have the following theorem. 

Lemma 2 (Topological shrinkage). For D- dimensional STS models (D — 1, 2, 3), let R and R' be reference regions: 
R,R' e Rref- When R' - R, one can deform geometric shapes of logical operators continuously from R' to R: 

ga^g^ for R^W. (77) 

Proof. The proof of the lemma is straightforward from theorem [T] and theorem |5] so we present a proof only for 
Ri - R\. Since all the zero-dimensional and one-dimensional logical operators can be supported inside we have 
gRj > k. Similarly, we have > k. Then, since ^r, + g-^ - 2k, we have ^r, = g-^ - k. □ 

Topological deformations of logical operators: So far, our discussion on a topological property of logical oper- 
ators has been limited only to regions generated by m-dimensional unit regions while some connected regions (such 
as examples presented in Fig. [24] ) cannot be generated by taking unions of unit regions. Then, a naturally arising 
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question is whether one can perform similar continuous deformations for arbitrary connected regions of composite 
particles or not. Here, we provide a complete description on topological properties of logical operators in STS models 
by extending the notion of continuous deformations to arbitrary connected regions. 

It turns out that one can continuously deform a geometric shapes of a logical operator defined on any connected re- 
gion of composite particles continuously for STS models with D = 1, 2, 3. Here, we summarize topological properties 
of logical operators in STS models as follow^ 

Theorem 6 (Continuous deformation). Consider STS models for D — 1,2,3. Consider two connected regions of 
composite particles R and R' which are topologically equivalent: R ^ R'. Then, we have 

gR = gR, for R - R'. (78) 

Therefore, for any given logical operator { defined inside R, one can always find an equivalent logical operator defined 
inside R' as long as R can be continuously deformed into R' . 
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Figure 24: Examples of continuous deformations. Note tliat systems have periodic boundary conditions. 



For clarity of presentation, we skip the proof of theorem |6] One may verify it by using scale symmetries, a bi- 
partition theorem (theorem [TJ, and the translation equivalence of logical operators (theorem [2|. We present some 
examples of continuous deformations in Fig.[24| 



6. Summary and open questions 

In this paper, we have established the connection between the feasibility of self-correcting quantum memory and 
the thermal stability of topological order, and provided partial answers to these two open problems by solving a model 
that may cover a large class of physically realizable quantum codes. Our discussion is limited to stabilizer codes with 
translation and scale symmetries, and these two problems still remain open for an even larger class of gapped spin 
systems on a lattice. Yet, we hope that our analysis will provide an important insight and a useful guidance on studies 



'Mathematically inclined readers may want to have more precise definitions of geometric shapes and continuous deformations. Here, we make 
some comments on this issue to make our discussion more rigorous. For a D-torus with ni X ■ ■ • X «£> composite particles, we split the entire system 
into ni X ■ ■ ■ X «o hypercubic regions where each hypercube contains each composite particle, and cover the entire system completely, but without 
any overlap. For a given set of composite particles R, we define its geometric shape as a union of hypercubes which contain all the composite 
particles inside R. After assigning geometric shapes to sets of composite particles in this way, one can define a continuous map between two sets R 
and R' , and infi'oduce an equivalence relation between them in a straightforward way. 
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of coding and physical properties of gapped spin systems on a lattice. We also hope that our work will contribute to 
stimulating the use of quantum coding theoretical concepts in studying many-body quantum systems further. 
Below, we discuss possible future problems and make some comments on them. 

Frustration-free systems with scale symmetries: While our discussion is limited to stabilizer codes, any gapped 
spin systems with degenerate ground states can be used as quantum memory devices in principle. It seems that, for an 
arbitrary gapped spin system defined on a lattice, there exists some frustration-free Hamiltonian which approximates 
the original system and serve as its low energy effective Hamiltonian (except some subtle properties such as chirality). 
In fact, this claim has been rigorously proven for one-dimensional gapped spin systems (5F\ by showing that any 
ground state of one-dimensional gapped spin Hamiltonians can be efficiently simulated through the matrix product 
state formalism. Now, this claim seems to be widely believed among the condensed matter physics community, and it 
is at the heart of recent progress in classifications of quantum phases in gapped spin systems lfT7ll57ll58l . Therefore, 
the analyses on coding properties of arbitrary frustration-free Hamiltonians with translation and scale symmetries may 
provide a useful insight on questions concerning coding properties and quantum phases arising in arbitrary gapped 
spin systems on a lattice. 

To the best of our knowledge, all the examples of frustration-free Hamiltonians with translation and scale symme- 
tries, such as the quantum double model lU and the string-net model 1591 , have the dimensional duaUty on geometric 
shapes of certain operators which may be considered as generalizations of logical operators. Also, these operators 
seem to have the continuous deformability in a way similar to logical operators in STS models. With these ob- 
servations, we feel that our results are universal for all the frustration-free Hamiltonians with translation and scale 
symmetries, and thus, effectively true for arbitrary gapped spin systems defined on a lattice with a small number of 
ground states. Also, these observations imply that any gapped spin systems with scale symmetries (or a small number 
of ground states) may be described through TQFT, as discussed in Appendix B Yet, the connection between TQFT 
and coding properties of gapped spin systems on a lattice with scale symmetries must be further established. 

Beyond scale symmetries: While our treatments in the present paper are limited to quantum codes with scale 
symmetries, there are several interesting models of quantum codes which do not have scale symmetries Il32l |60l - l62l . 
Of particular interests are models proposed in f6V\ and in f3T] which are classical and quantum memories respectively 
with partially self-correcting properties. These models do not have scale symmetries since the number of logical 
qubits kfi is highly sensitive to the system size n, and there is no constant upper bound on kfj. These models are known 
to have logarithmically large energy barrier AE ~ LOG(L) with a large number of energy local minima. As a result, 
these models seem to have POLY(L) relaxation time with slow dynamics, which may result in POLY(L) bit or qubit 
storage time. 

Whether these models with broken scale symmetries are useful as classical and quantum memories is a compli- 
cated problem. First of all, for a system to be a efficient quantum memory device, it is desirable to have exponentially 
long qubit storage time: t ~ EXP(L), since it needs at least d gate operations to write or readout a logical qubit, and 
these encoding and decoding processes takes at least polynomially long time t ~ POLY(L). Also, it seems difficult to 
find an efficient decoding algorithm with an efficient error-correction for these models. Finally, these models undergo 
phase transitions at T = 0, which may imply their potential thermal instability. At this moment, there are many is- 
sues to be analyzed concerning implementability of quantum codes with broken scales symmetries from engineering 
viewpoints. 

Beyond TQFT: Stabilizer codes with broken scale symmetries are remarkable examples which may be beyond 
the description of the standard TQFT. Such models do not have continuously deformable logical operators, and are 
not expected to be characterized by topological properties of logical operators as in STS models. This observation is 
consistent with the fact that systems described by TQFT, in a sense of the axiomatic formulation developed by Atiyah, 
are allowed to have only a finite number of degenerate ground states. Thus, finding an effective field theoretical 
description for stabilizer codes with broken scale symmetries may be an interesting future problem which may lead 
to discoveries of novel quantum phases that are currently unknown and are governed by some deeper mathematical 
formalism than topology. 

To the best of our knowledge, all the examples of stabilizer codes with broken scale symmetries still possess 
certain discrete scale symmetries. STS models have continuous scale symmetries since kjj - k for all the system sizes 
«, while stabilizer codes with broken scale symmetries seem to have discrete scale symmetries since the number of 
logical qubits kjf has good scaling properties under global scale transformations: (ni, «2, ■■■)—> {cn\,cn2, ■ ■ ■ ) where 
c is some appropriate integer with behaviors such as k^fi - ck^j. The distinction between continuous and discrete 
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scale symmetries becomes particularly important when one considers effects of RG transformations. As is pointed out 
in ll63l . systems with continuous scale symmetries correspond to fixed points of RG transformations, while systems 
with discrete scale symmetries correspond to limit cycles of RG transformations. Physical realizations of systems 
with discrete scale symmetries include the Effimov effect which is recently of particular interest in the ultracold atom 
physics community. This observation implies that stabilizer codes with discrete scale symmetries correspond to limit 
cycles of RG transformations. Yet, the connection between limit cycles and discrete scale symmetries for lattice 
systems must be further established since discussion in ll63l is given primary for continuum systems. 

Spin glass behaviors and translation symmetry breaking: Another interesting property of stabilizer codes with 
broken scale symmetries is a glassy behavior with slow relaxation dynamics. This glassy behavior may be understood 
as a direct consequence of broken scale symmetries, and associated broken translation symmetries in the ground space. 

In general, even if a Hamiltonian possesses translation symmetries, its ground states may not possess translation 
symmetries. Let us assume that a stabiUzer Hamiltonian is translation invariant under unit translations of composite 
particles: Tm{H) - H for all m. Then, ground states of this stabilizer Hamiltonian may break translation symmetries: 
TmiVI^)) + W)- Some examples of systems with broken translation symmetries are presented in fTT '38', '641, and 
translation symmetry breaking of ground states arising in topologically ordered systems are studied in [38, 64|. 

An interesting connection between scale symmetries and translation symmetries is that when a stabilizer code 
has scale symmetries: k^^ - k, all the ground states are invariant under unit translations of composite particles: 
TmiVI')) - Vjj) regardless of the system size, as proven in U7J . In other words, if a system is coarse-grained such that 
scale symmetries are satisfied, translation symmetries are protected inside the ground space. On the other hand, if a 
system does not have scale symmetries, translation symmetries may be broken inside the ground space. In particular, 
when the system has a large number of logical qubits, translation symmetries are strongly broken such that there is no 
finite translation which keep ground states equivalent. 

This strong breaking of translation symmetries seems to be the key to the glassy behavior observed in stabilizer 
codes with broken scale symmetries. In conventional spin glass models, a Hamiltonian consists of mutually com- 
muting terms whose signs are random, and translations symmetries of the Hamiltonian are initially broken due to the 
randomness. This randomness in the Hamiltonian gives rise to spin configurations of ground states which are not uni- 
form over real space, and leads to a complicated and slow relaxation dynamics in spin glasses. On the other hand, in 
stabilizer codes with broken scale symmetries, stabilizer Hamiltonians possess translation symmetries initially. How- 
ever, spin configurations are not uniform since translation symmetries are strongly broken inside the ground space due 
to broken scale symmetries. Thus, systems with broken scale symmetries are expected to exhibit spin glass behaviors 
with slow relaxation dynamics and a large number of local minima. Yet, the connection between the glassy behavior 
and broken scale symmetries must be further established since the observation given here is very heuristic. 



Appendix A. Topological order at finite temperature 



In the main part of the paper, we have discussed the thermal stability of topological order by using expectation 
values of logical operators. In this appendix, we verify the use of expectation values of logical operators as topological 
order parameters, and define the thermal stability of topological order more rigorously. 

In Appendix A. 1 we begin by verifying that expectation values of logical operators can be used for topological 
order parameter at T = 0. In [Appendix A.2| we define topological order at finite temperature by extending the 
definition of topological order at zero temperature to finite temperature. In Appendix A. 3 we discuss whether the 
existence of a large energy barrier is sufficient for the thermal stability of the system or not. 



Appendix A.l. Topological order parameters and logical operators 

While the ground state properties change only slightly under small perturbations in topologically ordered systems, 
the original ground state properties will be eventually lost under large perturbations. One may see whether topological 
order is lost by checking whether a new ground state can be approximated by the original ground state through a local 
unitary transformation. However, this naive approach will require a substantial amount of computations and cannot 
capture changes of coding and physical properties of the ground state under perturbations. 

Fortunately, the loss of the ground state properties under perturbations can be quantitatively characterized by 
some physical quantities, called topological order parameters. Here, we characterize the stability of topological order 
against perturbations through topological order parameters. 
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Topological order parameter: Topological order parameters are physical quantities with which one can judge 
if topological order is lost under perturbations or not. In this light, one may deduce the necessary properties of 
topological order parameters. By its definition, a topological order parameter must be a global function of the system 
which does not change under local unitary transformations. Also, a topological order parameter must change only 
slightly when topological order is protected, while it undergoes some non-analytic change when topological order is 
lost as a result of perturbations. 

There have been several proposals for such topological order parameters. The number of ground states serves as 
a topological order parameter since the ground state degeneracy is protected for topologically ordered systems at the 
thermodynamic limit [35] . Another interesting proposal is to use a certain entanglement measure, called topological 
entanglement entropy El?], which corresponds to a constant correction to the entanglement area law. Since topolog- 
ical entanglement entropy is a non-local quantity which involves a large number of spins, it does not change under 
local unitary transformations. 

Below, we analyze the stability of topological order thorough expectation values of logical operators in order to 
further build the connection between topological order and quantum codes. 

Expectation values of logical operators: Let us first consider two-dimensional Toric code and denote a pair of 
anti-commuting logical operators as £ and r where { extends in the 1 direction and r extends in the 2 direction. Then, 
one may use the expectation values of the summations of translations of logical operators as order parameters: 

U{{) = \^Yj^I{{\ U(r) = ^Yj^ti'-) (A.1) 

V -V 

where L is the Unear length of the system. Here, we took the summation of logical operators in a symmetric way. 

Now, let us discuss if topological order is stable under a perturbation V by considering the following perturbed 
Hamiltonian with an initial bias: 

Hfie) = HsTs- €Ui{) + V. (A.2) 

Note that, when V = 0, the small bias breaks the ground state degeneracy and 

{UiOh^o = 1. (A.3) 

Recall that, in evaluating the expectation value of U{£) in a presence of V, one must take the limit where e goes to 
zero after taking the limit where L goes to infinity. Since the ground space is separated from excited states by a finite 
energy gap, one may neglect the effect of excited states. So, let us examine the effect of eU{{) and V on the ground 
space. As a result of eU{(), the ground state space is split into two sectors with £ - \ and £ - -I where the energy 
splitting between two sectors is e. On the other hand, the energy splitting induced by a perturbation V is exponentially 
suppressed by a factor exp(-L/Lo)- Then, in considering the effect of V on a sector with £ - I, one may neglect the 
effect of a sector with f - - I when e is (9(1). Now, we take the limit of L — » oo. Then, the ground space corresponds 
to a sector with £ = 1, and we have 

m£)},^o > (A.4) 

even when e goes to zero. However, if V is large enough, the ground state properties will be lost and {U(£))f^o will 
be close to the value for the ground state of V. In a similar way, one may consider the following initial bias: 

Hrie) = HsTs - eUir) + V (A.5) 

and have 

<t/(r)),^o > 0. (A.6) 

Therefore, the stability of topological order can be characterized by expectation values of U{£) and U{r). 

This argument implies that the loss of the ground state properties may be captured through the stability of logical 
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qubits encoded in the ground space. The initial bias of eU (£) represents the initial encoding of a logical qubit. The 
expectation value of U{{) under a perturbation V represents how much information is protected in a presence of 
perturbations. 

Stability of topological order: Based on these observations, let us formulate the stability of topological order in 
terms of logical operators. Let us first generalize the definition of logical operators. We call operators U which satisfy 
the following conditions generalized logical operators: 

[Hq,U] = (A.7) 
U + e'"/ and |f/|max - 1 in the ground space (A. 8) 

where the action of U is defined inside the ground space of the original unperturbed Hamiltonian Hq. In other words, 
generalized logical operators are any operators which does not change the energy of the system, but acts non-trivially 
inside the ground space. Generalized logical operators are not c-numbers inside the ground space. For example, 
a scaled summation of translations of logical operators is a generaUzed logical operator, while projectors onto the 
ground space are not generalized logical operators. 

Based on generalized logical operators, we define the stability of topological order in the following way. 

Definition 3 (Stability of topological order). The system has topological order when there exists a pair of generalized 
logical operators IJ{1) and lJ(r) which satisfy the following conditions: 

• and C/(r) do not commute with each other: \U{C), t/(r)] + inside the ground space. This implies the 
existence of the ground state degeneracy, and U (€) and U (r) may characterize a logical qubit or qudit. 

• For perturbed Hamiltonians: 

He = Ho-eU({) + V, Hr = //„ - et/(r) + V (A.9) 
with any types of local perturbations V, 

{U({)),^o,v^Q = 1 for He (A. 10) 

<f/(r)),^o, v^o = 1 fori/,. (A. 11) 



In evaluating {U{{))e^o_v^Q, we first take the thermodynamic limit, and then, take the limit of V — > 0, and finally, 
take the limit of e — > 0. It should be emphasized that we consider any types of perturbations V in analyzing the 
stability of topological orde:[^ 

Let us look at some examples here. For a one-dimensional classical ferromagnet, let { and r be a zero-dimensional 
and a one-dimensional logical operators: 

^ = Zi, r = W^j- (^-12) 

i 

Then, for V - t YijZj, we have 

^ ^ (A. 13) 

j 

At the thermodynamic limit, we have 

{U({)) = -1 for He (A. 14) 



**The definition used iiere is similar to tlie one used in I44i except for taking tlie limit of V — > 0. This is because the definition in 1441 is 
appUcable to systems with gapless energy spectrum, while our definition is applicable only to systems with a finite energy gap. 
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regardless of f > and e > 0. Then, taking the Umit of f — > and e — > 0, we have 

{U{()),^o,v^Q = -1 foi-Hf. (A. 15) 

Therefore, the system is not topologically ordered. 

Topological order parameters and local unitary: Strictly speaking, expectation values of logical operators can- 
not be used as topological order parameters. In fact, while logical operators are non-locally defined for topologically 
ordered systems, their expectation values may change under local unitary transformations. However, topological order 
must be characterized by some physical quantities or objects which are not affected by local unitary transformations. 
In fact, both the ground state degeneracy and topological entanglement entropy do not change under local unitary 
transformations. In this light, one might think that the above definition of the stability of topological order is not 
legitimate. Here, we make some comments on this issue briefly. 

Despite the fact that expectation values of logical operators are not topological order parameters, we can formu- 
late the stability of topological order through them. This is because we considered not only a specific type of local 
perturbations but all the types of local perturbations. Here, we demonstrate that, if we consider only one type of per- 
turbations, expectation values of logical operators may fail to capture quantum phase transitions. In two-dimensional 
Toric code, let us consider a one-dimensional logical operator £ which consists only of Pauli X operators. Then, 
consider the following perturbation: 

y = (A. 16) 

where t is some positive parameter and r represents the position of each qubit. Then, when we increase f, there 
must be a quantum phase transition and the ground state properties will be close to the ones for V. However, such 
a perturbation may not change the value of {U({)) since V commutes with the logical operator £. In general, if one 
carefully choose the types of perturbations, expectation values of logical operators remain unchanged. This implies 
that if one tries to study a quantum phase transition only through {U(£)), one may fail to detect the transition]^ 

Appendix A.2. Topological order at finite temperature 

So far, we have focused on the ground state properties of topologically ordered systems at zero temperature 
and seen that topological order characterizes the ground state properties which are stable against any types of small 
perturbations. Here, we address the thermal stability of topological order. 

The stability of topological order against perturbations can be discussed by seeing if the perturbed ground state 
can be reached by local unitary transformations or not. However, the stability of topological order against thermal 
fluctuations cannot be formulated in a similar way since there is no unitary transformation which connects a pure 
state and a statistical ensemble. Here, we define the stability of topological order against thermal fluctuations by the 
changes of topological order parameters. 

Stability against thermal fluctuations: Below, we give the definition of the stability of topological order through 
expectation values of logical operators. 

Definition 4 (Topological order at finite temperature). The system is said to have topological order which is stable at 
finite temperature if and only if there exists a pair of non-commuting generalized logical operators U({) and U{r), and 
some finite transition temperature which satisfy the following conditions: 

• Consider perturbed Hamiltonians with initial biases 

H( = Ho-€U({) + V, Hr = Ho-£U(r) + V. (A. 17) 



'while expectation values of logical operators cannot be used as topological order parameters, geometric shapes of dressed logical opera- 
tors 1171 . which are operators characterizing the transformations inside the ground space of the perturbed Hamiltonians |55 |, may be used. Since 
continuous deformations between topologically distinct logical operators are not allowed, geometric shapes of dressed logical operators must 
undergo some discontinuous change during the transition. Thus, some non-analyticities will be induced in ground states during quantum phase 
transitions between Hamiltonians with topologically distinct logical operators. 
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For r = 0, 



(f/(r)),^o. v^o = 1 



for//f 
for Hr, 



(A. 18) 
(A. 19) 



and for Tc>T >Q, 



for//f 
for//. 



(A.20) 
(A.21) 



for any types of V. 

One may see the connection between self-correcting quantum memory and the stability of topological order. In 
particular, if topological order in a spin system is stable at finite temperature, such a system works as self-correcting 
quantum memory. A logical qubit encoded with respect to U{{) and U(r) can be read out by measuring U({) and U{r) 
even at finite temperature. This implies that encoded logical qubits are not lost in a presence of the interaction with 
the external environment, and there must be some self-correcting thermal dissipation processes. 

It should be noted that the expectation values of logical operators i and r usually vanishes at any finite temperature 
as one may see from a direct calculation. Only the symmetric summations of logical operators may give rise to non- 
vanishing expectation values of generalized logical operators. 

Comment on the definition of stability: Our definition of the stability of topological order relies on topological 
order parameters associated with a pair of non-commuting logical operators £ and r. The definition of the stability of 
topological order in the present paper is motivated purely from quantum information theoretical viewpoints. However, 
the definition we have used in the present paper may be too strict. For example, if one is interested only in the ground 
state properties associated to the expectation value of it makes perfect sense to claim that a three-dimensional STS 
model has topological order which is stable at finite temperature. In fact, once the system is held at finite temperature, 
the system properties are stable against any types of small perturbations. 

Similarly, the stability of the ground state properties depends crucially on types of perturbations. For example, 
while there is no topological order in a one-dimensional system in a strict sense, the ground state properties of a one- 
dimensional classical ferromagnet are stale against small perturbations if one declares to consider only perturbations 
which are products of X operators. Similarly, the ground state properties of the one-dimensional AKLT model is stable 
against small perturbations which do not break the time-reversal symmetry. Therefore, in discussing the stability of 
topological order, one needs to specify types of symmetries which are of interest. 

Appendix A. 3. Thermal stability of topological order and energy barrier 

We have seen that a large energy barrier in classical or quantum memory is the key to the thermal stability of 
ferromagnetic order or topological order at finite temperature. Yet, it is not clear if the existence of a large energy 
barrier is sufficient for the thermal stability or not. 

While STS models always have energy barrier 0{L") where a is an integer (a > 0) due to the dimensional duality of 
logical operators, there are several examples of stabilizer codes with broken scale symmetries which have logarithmic 
energy barrier: A/i ~ LOG(L). Let us look at an example of classical memory with broken scale symmetries [|6T1 . 
The model is constructed on a square lattice with L x L qubits: 



The model does not have scale symmetries, and there is no upper bound on the number of logical qubits (bits) k since 
k is highly sensitive to L. The model has logarithmic energy barrier AE ~ LOG(L), and as a result, the bit memory 
time is expected to be t ~ POLY(L) if one trusts the Arrhenius law. 

Despite a large energy barrier which scales with the system size, the model is known to be thermally unstable. This 
may be easily verified by computing the partition function: Z(J3) = Tre^^^ . Since N - k stabilizers are independent. 




(A.22) 
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we have 



(e-P + ^f-\e-Pf < Zifi) < {e-P + ^f-\^f (A.23) 

and, we have 

hm ^logZOS) = log(e"^ + /). (A.24) 

Therefore, the thermodynamic property of this model is equivalent to a single qubit in a magnetic field, and the 
model does not have the thermal stability. A similar discussion holds for stabilizer codes with qubits, supported 
by interaction terms, which include one-dimensional Ising model, two-dimensional Toric code, models proposed 
in ||321|60| . Note that M > in two-dimensional Ising model and four-dimensional Toric code where M is the number 
of interaction terms. 

It is not clear whether polynomial energy barrier is necessary for thermal stability of ferromagnetic order or not. 
Another important characteristic of models with logarithmic energy barrier is the existence of a logarithmically large 
number of local minima. As a result, at finite temperature, the entropic term easily dominates the free energy function. 
So, it seems that the thermal instability results from both the logarithmic energy barrier and the logarithmically large 
number of local minima. 

It is worth mentioning that there is an interacting spin model with a logarithmic energy barrier, but is not defined 
on a lattice. The model is called two-dimensional XY model, which is known to undergo a thermal phase transition at 
finite temperature. The model is of particular interest since it was not expected to possess any thermal phase transitions 
as a result of the Mermin-Wagner theorem which states that continuous symmetries cannot be spontaneously broken 
at finite temperature for D < 2. The reason why two-dimensional XY model undergoes a thermal phase transition 
at finite temperature is because the transition, known as Kosterlitz-Thouless transition, is induced by topological 
defects which are not directly related to continuous symmetries of the model. In two-dimensional XY model, it takes 
logarithmic energy to create a vortex: AE ~ log r where r is the size of a vortex. As a result, at low temperature, a 
configuration without vortices is favored, while at high temperature, a configuration with a large number of vortices 
is favored. Between these configurations, a phase transition occurs. 

From the observations above, it seems that the logarithmic energy barrier may lead to the thermal stability, while 
the existence of local minima may break the thermal stability. Yet, at this moment, the connection between the energy 
barrier and the thermal stability has not been completely established. 



Appendix B. Topology and quantum codes 

In section|5] we have shown that one can deform geometric shapes of logical operators continuously while keeping 
them equivalent. This continuous deformability of logical operators implies that one can introduce the notion of 
topology in analyzing and classifying coding and physical properties of STS models. In this appendix, we further 
discuss the role of topology in analyzing stabilizer codes. 

Many interesting properties of STS models, such as the dimensional duality of logical operators, naturally appear 
as corollaries of the continuous deformability of logical operators. In [Appendix B.l| we begin by showing that the 
dimensional duality of logical operators can be derived only by assuming the continuous deformability of logical oper- 
ators. In particular, we show that, for stabilizer codes with continuously deformable logical operators, m-dimensional 
and D - m-dimensional logical operators always form anti-commuting pairs where m is an arbitrary positive inte- 
ger (m < D). As an example of D-dimensional stabilizer codes with continuously deformable logical operators, we 
present generalizations of the Toric code to D-dimensional systems with anti-commuting pairs of m-dimensional and 
D - m-dimensional logical operators for arbitrary positive integers D and m. 

The topological property of logical operators naturally leads us to consider a possible relevance to another well- 
celebrated theoretical framework equipped with the notion of topology; called topological quantum field theory 



(TQFT) 116514671 . In Appendix B.2 we establish the connection between quantum codes with the continuous de- 
formability and TQFT further by demonstrating that the braiding of anyonic excitations in a D-dimensional stabilizer 
code is characterized by a topological invariant in a D + 1-dimensional system. This implies that such a D-dimensional 
stabilizer code can be effectively described by D + 1 -dimensional TQFT. 
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Appendix B.l. Dimensional duality as a corollary of continuous deformability 

The continuous deformability of logical operators implies that topology is the essential notion in analyzing coding 
and physical properties of STS models. Then, a naturally arising question concerns the role of topology in determining 
coding and physical properties of stabilizer codes. Here, we show that the dimensional duality of logical operators 
is a universal property for all the stabilizer codes with continuously deformable logical operators. We also give 
concrete examples of D-dimensional stabilizer codes which have continuous deformability of logical operators and 
the dimensional duality by generalizing the Toric code to D-dimensional systems. 

Dimensional duality from topological deformation: Let us begin by counting the number of independent logical 
operators defined inside m-dimensional regions. Recall that m-dimensional concatenated unit regions are obtained by 
taking unions of all the m-dimensional topological unit regions. Then, one may call logical operators which can be 
defined inside /?,„, but cannot be defined inside Rm-\, m-dimensional logical operators. 

Definition 5. m-dimensional logical operators are logical operators which have representations supported inside R,„, 
but do not have representations supported inside Rm-\- 

Now, let us denote the number of independent m-dimensional logical operators as g,„ = g^^^^ - g^^ ^ , where = guo 
by setting gR ^ = 0. Then, there exists an interesting relation among the numbers of m-dimensional logical operators. 
In particular, the following lemma holds. 

Lemma 3. There are the same number of m-dimensional and D — m-dimensional logical operators: 

gm = gD-m for m = 0, ■ ■ ■ , D. (B.l) 

The proof of this lemma can be obtained through a simple algebra by combining theorem[T|and the deformability 
of logical operators. 

Proof. Consider a bi -partition of the entire system into R,„ and R„. From the topological deformation of logical 
operators, we have 

gR;: = gRD-„-i (B-2) 

since R„ ^ Ro-m-i- Thus, R„, and R,,, support the following logical operators: 

Rm : 0-dim, 1-dim, ■ ■ ■ , m-dim 

Rm : 0-dim, 1-dim, ■ ■ ■ , D - m - 1-dim. 



Therefore, we have 

m D-m- 1 

Recall that gR gR = 2k as presented in theorem[T| Using this formula for R - R,„, we have 

gR„^+g^^2k. (B.4) 

Then, we have 

m D-m- 1 

Sj^^k^ for m = 0,---,D. (B.5) 

M) .1=0 

Since the total number of independent logical operators is E)Lo gj - we have g„j - go-m for all m. □ 
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This lemma implies the existence of a dimensional duality in geometric shapes of logical operators. One may 
notice that this lemma is just a manifestation of Poincare duality in a D-torus where the mth and D - mth Betti 
numbers are equal. 

To completely establish relations between each logical operator with different dimensions, let us analyze their 
commutation relations. We have the following theorem. 

Theorem 7 (Dimensional duality). One can choose a set of 2k independent logical operators of D-dimensional 
systems with the deformability of logical operators in the following way: 



4 



(B.6) 



where tp are nip- dimensional logical operators and rp are D - mp- dimensional logical operators for some integer mp 
(Q < mp < D) for any p — I, - ■ ■ ,k. 

In other words, one can choose logical operators such that the summation of dimensions of pairs of anti-commuting 
logical operators is always D. Theorem |7] follows immediately from the following lemma. 

Lemma 4. m-dimensional and m' -dimensional logical operators commute with each other ifm + m' < D. 

Proof. Consider a m-dimensional logical operator { and a ;7i'-dimensional logical operator {' which is defined inside 
R,n and R„f respectively. For m-hm' < D, there exists a translation of R,„ such that R,„ and R„f have no overlap. Then, 
due to the translation equivalence of logical operators, some translation of £ do not have an overlap with t, which 
leads to [{, £'] =0. □ 

With this lemma, the proof of theorem |7] is immediate by using lemma [3] For example, from the lemma, zero- 
dimensional logical operators may anti-commute only with D-dimensional logical operators. Since there are the same 
number of zero-dimensional and D-dimensional logical operators, there exists a canonical set of logical operators 
where D-dimensional logical operators can anti-commutes only with zero-dimensional logical operators. Similarly, 
one can show that there exists a set of 2k independent logical operators such that m-dimensional logical operators 
anti-commute only with D - m-dimensional logical operators for all m. 

Generalization of the Toric code: Now, we give concrete examples of D-dimensional stabilizer codes which 
have continuously deformable logical operators. The model we present here is a straightforward generalization of 
two-dimensional Toric code to D-dimensional settings for arbitrary positive integer D. In particular, we illustrate the 
construction of D-dimensional Toric code with anti-commuting pairs of m-dimensional and D — m dimensional logical 
operators for arbitrary positive integers D and m (m < D). 



Star 




plaquette 




Figure B.25: An example of the construction for D = 3 and m = 2. 



We first consider a D-dimensional hypercubic lattice which consists of = L x ■ ■ ■ x L of D-dimensional unit 
hypercubes with periodic boundary conditions. We denote a set of /7-dimensional unit hypercubes in this lattice as 
hp. Note that \hp\ - oCp ■ N since one needs to specify p directions from D directions in defining /^-dimensional unit 
hypercubes. We put qubits at the centers of m-dimensional hypercubes. Then, the total number of qubits is oCm ■ N, 
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and qubits are labeled by m-dimensional unit hypercubes in h,„. Fig. B.25 shows a construction of the model for D = 3 
and m-2 where qubits live at centers of two-dimensional unit squares. 

The Hamiltonian consists of plaquette terms and star terms as in the conventional two-dimensional Toric code: 

^ ~ ^plaquette H star 

(B.7) 

Here, we call (m -i- l)-dimensional unit hypercubes in h,„+i "plaquettes". Then, plaquette terms consist of Pauli Z 
operators which act on qubits included in plaquettes: 

Hplaquelle - — ^ Bp, Bp - ]"[ Zr (B.8) 

peh„,+i rep 

where r represent m-dimensional unit hypercubes included inside a m + 1 -dimensional unit hypercube p. Note that 
there are oCm+x ■ N plaquette terms. Next, we call (m - l)-dimensional unit hypercubes in hm-i "stars". Then, star 
terms consist of Pauli X operators which act on qubits neighboring to a star: 

Hstar = - 2 A„ A, = (B.9) 

where s c r means that a star s is included inside m-dimensional unit hypercube r. Note that there are oCm-i ■ N star 
terms. Noting that [As, Bp] = since As and Bp share either zero or two qubits in common, the model is a stabilizer 



code. Fig. B.25 shows constructions of star terms and plaquette terms for D - 3 and m-2. 

One can verify that the model has k - dC„, logical qubits with k anti-commuting pairs of m-dimensional and (D - 
m)-dimensional logical operators where a m-dimensional logical operator consists of Pauli Z operators supported on a 
m-dimensional hyperplane, while a (D - m)-dimensional logical operator consists of Pauli X operators supported on a 
(D - m)-dimensional hyperplane. Since they may share either zero or one qubit, they may commute or anti-commute 
with each other. One may easily see that the construction above reproduces two-dimensional, three-dimensional and 
four-dimensional Toric code for choices of {D, m) = (2, 1), (3, 1), (4, 2). For m - and m = D, the model is reduced 
to the D-dimensional Ising model. One can verify that logical operators arising in this model can be deformed 
continuously by using the bi-partition theorem (theorem [TJ. The construction above can be easily generalized to 
arbitrary D-dimensional graph embedded in a D-dimensional geometric manifold in a way similar to two-dimensional 
Toric code. 



Appendix B.2. D-dimensional STS model and D + l-dimensional TQFT 

The emergence of the notion of topology in geometric shapes of logical operators leads us to consider a possible 
relevance to topological quantum field theory (TQFT), which also deals with systems whose physical properties 
depend heavily on topological characteristics of the systems. Here, we make an attempt to establish the connection 
between stabilizer codes with continuously deformable logical operators and TQFT further 

Roughly speaking, TQFT is a field theory which is invariant under diffeomorphism (continuous deformation), 
and particularly suited for describing topologically ordered systems. The most important characteristic in systems 
described by TQFT is the invariance of all the correlation functions under diffeomorphism. Consider an arbitrary 
diffeomorphism through a continuous change of space-time coordinates x — > x'. Then, the correlation function of a 
scalar operator <p(x) satisfies 

{0i\</>(xi)cl>(x2) ■ ■ ■ <f>(x„)\Oj) = <0#(x'i),^(4) ■ ■ ■ <p(x'„)\Oj) (B. 10) 

where |0,) represent degenerate ground states. Therefore, the vacuum expectation value of any products of scalar 
operators is invariant under differmorphism, and only the topological properties of products may characterize their 
expectation values. 

Below, in order to establish the connection between stabilizer codes and TQFT, we show that the braiding of 
anyonic excitations in D-dimensional stabilizer codes can be characterized by some topological invariant in a D -H 1- 
dimensional system. In particular, by characterizing propagations of anyonic excitations in a D + l-dimensional 
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system, we demonstrate that the braiding of anyonic excitations corresponds to a configuration of m-dimensional and 
D - m-dimensional closed objects with a non-zero linking number 

We begin by analyzing propagations of anyonic excitations in two-dimensional Toric code. Anyonic excitations 
can be created by applying a segment of a one-dimensional logical operator to a ground state of the Toric code 
Hamiltonian since endpoints of a segment of a logical operator f"^ may not commute with interaction terms and 



create localized excitations (Fig. B.26 a)): i'"'^\i//gs), and one may make anyons propagate along a geometric shape of 



a one-dimensional logical operator. Since one can deform a geometric shape of a one-dimensional logical operator 
continuously in the Toric code, anyons can propagate freely on the lattice by applying a segment of a deformed one- 



dimensional logical operator (Fig B.26 b)). Therefore, the continuous deformability of logical operators is the key to 
propagations of anyonic excitations. 



(a) 



(b) 



anyons 



anyons' 



Figure B.26: Anyonic excitations created from segments of logical operators, (a) A one-dimensional logical operator, (b) A deformed one- 
dimensional logical operator. 

Propagations of anyonic excitations can be characterized by a one-dimensional closed loop drawn in a three- 
dimensional space. To illustrate this point, let us consider a process of creation, propagation, and annihilation of 



anyons, as described in Fig. B.27 a). By drawing propagations of anyonic excitations in a three-dimensional system by 



adding the time axis, the entire process can be represented as a one-dimensional closed loop as shown in Fig. B.27 b). 



In general, arbitrary one-dimensional closed loop in a three-dimensional system may characterize some propagations 



of anyonic excitations. See an example in Fig. B.28 which involves creations of two pairs of anyonic excitations, 
(a) creation propagation 



annilnilation 







Figure B.27: The correspondence between anyonic excitations in two-dimensional Toric code and a closed loop in a three-dimensional system, (a) 
Creation, propagation and annihilation of anyonic excitations, (b) Anyonic excitations described in a three-dimensional system. 

Next, let us consider the braiding between anyonic excitations in the Toric code. One interesting property of 
topologically ordered spin systems is the non-trivial braiding property between anyonic excitations. There are two 
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creation 



Figure B.28: Propagations of anyonic excitations and an associated closed loop in a three-dimensional system. 



distinct type of anyons in two-dimensional Toric code which are created by a pair of anti-commuting logical operators 
respectively, and the braiding of different types of anyonic excitations may give rise to non-trivial change inside the 
ground space. In two-dimensional Toric code, the braiding of different anyons give rise to an additional phase -1 
to the original ground state li/'j.J — > -\i//gs), which is a direct consequence of the anti-commutation between logical 
operators. 

One may understand this non-trivial braiding arising in two-dimensional Toric code as a topological invariant in 
a three-dimensional system. Let us consider the braiding of anyons described in Fig. B.29 a) where two types of 



pairs of anyonic excitations are created. One can characterize this braiding process as two closed loops in a three- 



dimensional system where loops are linked as described in Fig. B.29 b). In a more technical language, the braiding 
of anyons occurs only if the linking number between two loops has non-zero value. In particular, the final state is 
\^gs) (-l)'^'"'' |iAg.s) where Nii„k is the linking number of a given configuration of one-dimensional closed loops. This 
observation indicates that the braiding of anyonic excitations can be characterized by topological invariants, such as 
the linking number, in a (2 -H l)-dimensional system. 



(a) 



o 



o 



o 



o 



o 



oo 

i 



(b) 




Figure B.29: The braiding as a topological invariant in a three-dimensional system, (a) A braiding of anyonic excitations, (b) Loops with non-zero 
linking number. 



A similar observation holds for D-dimensional Toric code with a pair of m-dimensional and D - m-dimensional 
logical operators. The first type of anyonic excitations can be created by a segment of m-dimensional logical operator, 
while the second type of anyonic excitations can be created by a segment of D - m-dimensional logical operator. 
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The propagation of anyonic excitations of the first type can be characterized by a m-dimensional object in a D + 1- 
dimensional systems, while the propagation of anyonic excitations of the second type can be characterized by a 
D - m-dimensional object. These anyonic excitations are braided when the linking number between m-dimensional 
and D - m-dimensional objects is non-zero. While we have discussed the braiding for D-dimensional Toric code, 
similar discussion holds for any stabilizer codes with continuously deformable logical operators. Therefore, one 
may expect that D-dimensional stabilizer codes with continuous deformability can be effectively described hy D + 1- 
dimensional TQFT. 



Appendix C. Decomposition of logical operators 

We present the proof of theorem [5] in this and the next appendices. The goal of this appendix is to prove the 
following theorem which will be the key to the proof of theorem |5] 

Theorem 8 (Decomposition). Consider a three-dimensional STS model with the system size ni — 2 ■ 2^"^''!, «2 — 2™ 
and arbitrary nj, where m is an arbitrary positive integer For a given logical operator I supported inside P(n\ , «2, 1). 
one can decompose t as a product of the following centralizer operators 



where 



and {a is defined inside P{2v, n2,l)- 



T':(€h) = 4, P < 2^"^" (C.2) 



Here, Cr represents the restriction of the centralizer group C onto a region of composite particles R, meaning 
that Cr is a subgroup of centralizer operators defined inside R. Therefore, ^,,4 e Cp(„,,„,,i) means that and ^ 



are centralizer operators defined inside P{n\,n2, 1). We show the claim of the theorem graphically in Fig. C.30 The 
theorem claims that a two-dimensional logical operator defined inside P(n\,n2, 1) can be decomposed as a product 
of a one-dimensional centralizer operator and a two-dimensional centralizer operator (h which is periodic in the 1 
direction. As we shall see later, "2v" comes from the number of independent generators for the Pauli group acting on 
a single composite particle. 

Before starting the proof of theorem |8] let us describe the entire sketch of the proof of theorem |5] As a simple 
extension of theorem [8] one can show that a one-dimensional logical operator defined inside P(2v, «2, 1) can be 
further decomposed as a product of a one-dimensional and a zero-dimensional centralizer operators. After these 
decompositions, one can classify geometric shapes of logical operators according to their dimensions and can find 
commutation relations between them. 

Although theorem [s] is Umited to some specially chosen system sizes: n\ - 2 ■ 2^"-' ! and = 2"', one can 
construct logical operators for arbitrary system sizes from theorem [8] For example, due to scale symmetries, one 
can show that one-dimensional logical operators found in theorem [Sjare also logical operators for the systems with 
arbitrary n\. In fact, one can find logical operators in the forms described in theorem |5] These arguments will be 



presented in Appendix D 



Appendix C.l. Sketch of proof of theorem^ 

First, we note that theorem [8] was proven for n2 - \ {m - 0) in fTT'l since such a system with n2 - \ can be 
considered as a two-dimensional system which extends only in the 1 and 3 directions. For a two-dimensional STS 
model, we have the following lemma. 

Lemma 5. Consider a two-dimensional STS model where n\ —2 - 2^''! and arbitrary n2- For a given logical operator 
€ supported inside P{n\, 1), one can decompose { as a product of the following centralizer operators 

t ~ iatb, ia,tb ^ Cp(„,,l) (C.3) 
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P(ni,n2, 1) 




periodic in 1 direction 



Figure C.30: The claim of theorem[8] One can decompose a two-dimensional logical operator as a product of a one-dimensional centralizer operator 
Ca and a two-dimensional centralizer operator {),. 



where 



/5 < 2 



2v 



(C.4) 



and ta is defined inside P(2v, 1). 



We present the claim of the lemma graphically in Fig. C.31 The lemma claims that a one-dimensional logical 
operator defined inside P{ni, 1) can be decomposed as a product of a zero-dimensional centralizer operator and a 
one-dimensional centralizer operator ii, which is periodic. 



(a) 



-F(ni,l) 



(b) 



periodic in 1 direction 



2v 



Figure C.31: The claim of lemma|5] One can decompose a one-dimensional logical operator as a product of a zero-dimensional centralizer operator 
and a one-dimensional centralizer operator (t,. 

A three-dimensional STS model may be viewed as a two-dimensional system if one considers 1 x n2 x 1 composite 



particles as a single composite particle which consists of vn2 qubits (see Fig. C.32 1. In other words, we view the entire 
system as a two-dimensional lattice of one-dimensional tubes. Then, as a direct consequence of the lemma above, we 
notice the following corollary. 

Corollary 1. A logical operator { considered in theorem^can be decomposed as a product of the following centralizer 
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operators 



/'(«!, "2,1) 



(C.5) 



where 



Tl{(b) = 4, where (3 < 2^"^" 



(C.6) 



and £a is defined inside P(2vn2, n2, 1). 



We present the claim of the corollary graphically in Fig. C.32 The corollary claims that a two-dimensional logical 
operator defined inside P(ni,n2, 1) can be decomposed as a product of a one-dimensional logical operator and a 
two-dimensional logical operator ^ which is periodic in the 1 direction. 

However, a one-dimensional logical operator described in corollary |2] is not one-dimensional in a strict sense 
since it is defined inside P(2"' ■ 2v, 2"', 1), and its "width" 2"' ■ 2v grows as m increases. On the other hand, fa described 
in theorem[8]is truly one-dimensional since its width is at most 2v. Therefore, we need to show that a logical operator 
defined inside P(2™ ■ 2v, 2'", 1) have an equivalent logical operator defined inside P(2v, 2"\ 1). 




periodic in 1 direction 



2™ ■ 2u 



Figure C.32: The claim of corollary[2] Tlie width of a one-dimensional logical operator increases as m increases. 

The rest of this appendix is dedicated to the proof of the following lemma. 

Lemma 6 (Shrinkage). For system sizes considered in theorem [S] a logical operator operator t defined inside 
P(x, 2™, 1) always has an equivalent logical operator £' which is defined inside P(x — 1,2'", 1) when 2v < x < n^. 

By using this lemma, one can shorten the width of ia from 2"' ■ 2v to 2v. 

Appendix C.2. Identity generating matrix 

Here, we discuss how to shorten the width of {„. In particular, we introduce a certain binary matrix which is 
essential in reducing the width of £„. 

Shrinkage in two dimensions: To give an intuition on how to shorten the width of let us first consider the case 
where m = 0. (So, this is a two-dimensional system with n2 - I, and instead of the "width", we use the "length"). As 
an example, consider the case when m - 0, x = 4 and ( is given by 

( = [a, B, C, AB] (C.7) 
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where { is defined inside P(4, 1, 1), and A, B and C are some Pauli operators. Then, consider the following logical 
operator: 



[A. 



B, AC, I, BC, ABC, 



ab]. 



(C.8) 
(C.9) 



Note that £" ~ { since E" is a product of three logical operators which are equivalent to each other due to the translation 
equivalence of logical operators. Then, we notice that following operators are centralizer operators: 



A = [a, b, ac\, ti = [bc. 



ABC, 



ab\ 



(C.IO) 



where t" - {\T^({2) since stabilizers in STS models are defined inside 2x2 composite particles and cannot overlap 
with ii and T^({2) simultaneously. Now, due to the translation equivalence of logical operators, we have 



f = ^1^2 = [a, B, AC] X [bc, ABC, ab] = [ABC, AC, BC]. 



(C.ll) 



Thus, a logical operator (, with the length 4, is shrunk into an equivalent logical operator {', with the length 3. 

An important observation is that one can form an identity operator / by taking a product of Pauli operators in f. In 
general, if the length jc of ^ is larger than 2v, one can always form an identity operator / by taking a product of some 
Pauli operators in { since there are 2v independent generators for single Pauli operators acting on a single composite 
particle. Now, let us consider { represented as 



= [Ui, U2, 



which is defined inside P(x, 1,1) where x > 2v. 
(0, • • • ,0) which satisfies the following condition: 



Then, there always exists a binary vector B - (Bi 



Now, let us take the following product of translations of £: 



(C.12) 

,B,) It 

(C.13) 



(C.14) 



Then, one may readily know that the xth entry of the above operator is /. From this operator, one can find two 
centralizer operators. By using them, one can readily shrink the length of f from x to x - 1. This trick is the key 
to the proof of lemma |5] Although the argument above works only when B has an odd number of 1 entries, a slight 
modification makes the shrinkage of £ possible when B has an even number of 1 entries. 

Identity generating matrix: Now, we consider more general cases with m > 0. Let us represent a logical operator 
{ defined inside a region P{x, 2™, 1) as a x x 2"' matrix whose entries are single Pauli operators: 



(C.15) 



where each Pauli operator Ujj acts on each composite particle. We also represent each column of { as follows: 



Ui 



(C.16) 
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Here, we denote a group of Pauli operators supported by a single column P(l,2"', 1) as !P™, and call it the column 
operator group and its elements column operators. Note that Uj e !P™ and !P™ ^ has 2™ ■ 2v independent generators: 
G(!P'«,) = 2" ■ 2v. 

When m - and x > 2v, we found a binary vector B with x components which characterizes how to form an 
identity operator from £. When m > and x > 2v, we can find an x x 2™ binary matrix B which characterizes how to 
form an identity operator from £. Here, we introduce the identity generating matrices as follows. 

Definition 6. Consider a logical operator { defined inside P{x, 2™, 1). 

• For a X X 2"' binary matrix B 



B 



we define the following operations: 



B 



1,2'" 



Bij = 0,1, 



(C.17) 



m = 



ij 

X 2'" 



(C.18) 
(C.19) 



'=1 >=i 

where ({B) is a product of translations of I taken according to B while ({B)x is the xth column of ({B). 
We call a binary matrix Z? identity generating matrix if and only if 

ro, ■ ■ ■ 01 



I and B 



0, 







(C.20) 



We assign parities to each column of a binary matrix B as follows: 

Par{B)i = Y^Bij (mod 2) (C.21) 

j 

where i = 1, ■ ■ ■ , x. We call a binary matrix B odd if and only if 

3/ s.t PariB)i = 1. (C.22) 



Therefore, when we form an identity operator, we considered translations of { both in the 1 and 2 directions. The 
identity generating matrix is said to be odd when there exists a column with an odd parity. 

Slirinkage tlirougli odd matrices: Note that there always exists an identity generating matrix when x > 2v. 
Then, with the existence of identity generating matrices, one might hope that the width of can be reduced until it 
becomes 2v in a way similar to the cases where m = 0. However, there is a caveat. In fact, only identity generating 
matrices with some special properties can be used for shrinking. In particular, we have the following lemma. 

Lemma 7 (Slirinkage tlirougli odd matrices). If there exists an odd identity generating matrix for £ defined inside 
P(x, 2"\ I), £ has an equivalent logical operator £' defined inside P(x — 1,2™, 1). 

Therefore, if there exists an odd identity generating matrix for any x with ni > x > 2v, one can complete the proof 
of lemma |6] Below, we present the proof of lemma |7] The existence of an odd identity generating matrix will be 
proven later 
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Proof. Assume that B is an odd identity generating matrix for {. Assume that for some /' (1 < /' < x), we have 



Par(B)i, = 1 and Par(B)i = for / < /'. (C.23) 
So, /' is the smallest integer such that /'th column has an odd parity. Here, we define the following binary matrix B' 



(see Fig. C.33 1 



= 



( i < i' ) 
( / > /' ). 



(C.24) 
(C.25) 



Note that B' consists of /th columns of B with / < /'. Based on B' , we consider the following logical operator (': 

£' = €{B') ~ t. (C.26) 

Note that €' is equivalent to { since {' is a product of an odd number of translations of (. (Note that {{B) may not be 
equivalent to I since the number of 1 entries in B may be even). 



See Fig C.33 for graphical representations of B, B', C{B) and C{B'). Note that ({B) has an identity operator at xth 
column. Note that ({B') has identity operators in the first x - /' columns since Bi j - for / > /'. 



even -i ,odd 




—r 

identity 



Figure C.33: Constructions of B, B' , C(B) and i(B'). 



Since t(B) has an identity operator at xth column, we can decompose it as a product of two centralizer operators 
whose lengths are at most x - 1. Let us denote these centralizer operators as JJa and f//,: 



{C.21) 



where Ua is the centraUzer on the left hand side and Ub is the centralizer on the right hand side (See Fig. C.33 1. Ua is 
defined from 1st column to x - 1th column, and Ub is defined from x + 1th column to 2x - 1th column. 

Since B and B' have the same entries from 1st column to /'th column, we notice that ({B') and Uh have the same 
Pauli operators from 2x - 2 - /'st column to 2x - 1th column as shown in Fig. C.33 Then, by applying Ub to ({B'), 

Although Ubt{B') may not be 



one can shrink the size of £{B'). In particular, Ubf{B') has the length at most x - 1. 
equivalent to {{B') as Ub may not be a stabilizer, one may consider the following logical operator: 

Ub X {(B') X T-''(Ub) ~ e(B') ~ t 



(C.28) 



which is equivalent to t{B') due to the translation equivalence of logical operators, and is defined inside a region with 
X - 1 X 2'" X 1 composite particles. Then, due to the translation equivalence of logical operators, there exists a logical 
operator t" ~ t which is defined inside P{x - 1,2™, 1). This completes the proof. □ 
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Appendix C.3. Existence of an odd matrix for m — I 

Next, we present a proof of the existence of an odd identity generating matrix for x > 2v, in order to complete the 
proof of lemma |6] and theorem [8] In particular, we shall prove the following lemma. 

Lemma 8. When x > 2v, there always exist an odd identity generating matrix. 

We start by discussing cases with m - and m - 1 before presenting general discussion. First of all, when m - 0, 
identity generating matrices are always odd since all the binary matrices are odd except B - (Q, - ■ ■ ,0). Therefore, we 
consider the cases where m - \ below. 

Characteristic value: Recall that we represented ^ as a jc x 2 binary matrix: 



f/l,2, 



and each column of I as follows: 



Uj,2 



U,,2 



U = I,-- - ,x). 



(C.29) 



(C.30) 



Now, we suppose that there is no odd identity generating matrix for f, in order to use the contradiction for the proof 
of lemma [8] 

First, it is convenient to classify column operators in into two types as follows. For a column operator U 
represented as 



U = 



Ui,2 



(C.31) 



we assign a characteristic value b and characteristic operator V as follows: 

• If = t/i,2, we assign a characteristic value b = I and a characteristic operator V - t/ij. 

• If Ui^i + Ui^2, we assign a characteristic value b = and a characteristic operator V - t/i,i f/i,2- 

One may easily understand this classification by representing U explicitly. A column operator U with = 1 is 



and a column operator U with b - Qis 



U 



U = 



VUi,2 
Ui,2 



(C.32) 



(C.33) 



So, a column operator with b - I is symmetric while a column operator with b - is not. Note that a characteristic 
operator is not an identity operator / except when U - I. 
Next, we introduce the following x x 2 binary matrices: 



£(i';0) = 



El.2, 



E,,i 

E,,2 



(C.34) 



such that 



Ei.i 
Eij 



a = i') 

otherwise 



(C.35) 
(C.36) 
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and 



£(/';!) = 



El.2, 



E,,2 



such that 



For example, 



and 



Eij = 



£(2;0) 



E{2-\) = 



otherwise 



0, 1, 0, 
0, 0, 0, 



0, 1, 0, 
0, 1, 0, 



(C.37) 



(C.38) 
(C.39) 



(C.40) 



(C.41) 



Now, let us represent characteristic values and characteristic operators for each column of £ as follows: 



bu Vi 

b2, V2 

b,, V,. 



(C.42) 



Then one can establish a connection between binary matrices E(i; 0) and E(i; 1), and characteristic values bj and 
operators Vj as follows: 



({E(i;0)), = 



when bi - 1 
when bi - 0. 



(C.43) 
(C.44) 



Proof of lemma |8f Now, let us proceed to the proof of lemma |8] for m - \. Without loss of generality, we can 
assume that 



b, = 1 



(C.45) 



for xq < X since permutations of columns do not affect the parities of identity generating matrices. (We will justify 
this later). We define the following sets of integers: 



b(0) = {I,-- - ,xo), b(l) ={x^, + \,■■■ ,x]. 

We denote groups of Pauli operators generated by V, with bj - and Vj with bj = 1 as 'Vo and 'Vi : 

^0 = < { Vi : i e b(0) ) > 
^1 = < { y,. : / e b(l) ) > . 



(C.46) 



(C.47) 
(C.48) 



Let us show that a set of characteristic operators {V,) for / e b(l) is independent : G{'V\) - x — xq. For this 
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purpose, we suppose that there exists some set of integers A c b(l) such that 



Then, the following binary matrix is an identity generating matrix: 

B = Y^E{i-0) 



(C.49) 



(C.50) 



/eA 







I 


Vi 




I 



(C.51) 



However, since B is odd with PariB), - 1 for / e A, this leads to a contradiction. 

Next, let us analyze 'Vq. For simplicity of discussion, first assume that [Vi] for i e b{Q) are independent. We 
consider more general cases where {V,) for / e b(0) are over complete later. Then, we have G{'Vq) - xq. Now, we 
define the following operators for / e b(0): 



t/; = UiT2m 
Vi ■ 

Note that U'. has a characteristic value b'. = 1 and a characteristic operator V,. Notice that 

f/; = e(E(i; 1)),. 

Since x > 2v, there exists a set of integer A such that 

]~[V,- = /, A ^ b(l) and A ^ b(0). 



(C.52) 



(C.53) 



(C.54) 



Note that A includes integers both from b(0) and b(l). Then, one notices that the following matrix B is an identity 
generating matrix: 



B = ^ £(/;!)+ ^ £(/;0) 



(C.55) 



/eAnb(O) 



leAnb(l) 



smce 



mx = n n 

(ieAnb(O)) |ieAnb(l)) 

= n 

igA 



V 




I 


Vi 




I 



(C.56) 



However, for an integer / in A n b(0), Par(B)i - 1, and B is an odd matrix. This leads to a contradiction. Note that 
discussion above is valid under the permutations of columns. 

Proof of lemma|8j continued: Next, let us consider the case where {V,) for / € b(0) are not independent. Let us 
denote the number of generators for 'Vo as io - GCVo) (xq < xq). Without loss of generality, we may assume that 
V\, - ■ ■ , Vf„ are independent since permutations do not affect parities of identity generating matrices. Here, we change 
our notations slightly: 



b(0) = {I,-- - ,.fo). 



b(0)' = {xo + 1,--- ,.xo), b(l) ={xo + l,--- 



(C.57) 
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and 

^0 = ({Vi : i € b(0) } > G(%) = xo (C.58) 
= {{Vi : i € h(l)}) Gm) = x-xo. (C.59) 

Since a set {V,} for bi = is over complete, there are xo - xq sets of integers A, (; = xq + 1, • • • , xq) such that 

Y\Vi' = I where i' < i, Si' e A, and i e A; (C.60) 

i'€Ai 

where the largest integer in A, is i. For i — xo + I, - ■ ■ , xq, we form the following operator: 



(C.61) 

which has a characteristic value b'- = I and a characteristic operator V^. Here, we denote a group of {V- } for j e b(0)' 



as 



% = ({y; : J G b(0)'}). (C.62) 
Now, we show that {V[} for i e b(0)' are independent: GCVq) = xq - xq. If there exists A c b(0)' such that 

f[ y; = /, (C.63) 



we have the following identity generating matrix: 



B = 2I]^0';0) (mod 2) (C.64) 



smce 



(€A 







7 






/ 



(C.65) 



Recall that the largest integer in A, is i. Let the largest integer in A be imax- Then, we have Par{B)i^ = 1, and thus, B 
is odd. This leads to a contradiction. 
So far, we have shown that 

GCVo) = xo, G(%) = xo-xo, GCVi) = x-xi. (C.66) 

Since x > 2v, there exists a set of integers A such that 

n n n = ^- ^^-^^^ 

{!€Anb(l)) {ieAnb(l)') {!€Anb(0)) 
The following matrix is the identity generating matrix: 

B = Yj Yj ^£0';0)+ Yj (mod 2) (C.68) 

i€Anb(0) !€Anb(0)' j'eA; !€Anb(l) 
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since 



n I n 



|/eA:i£b(0)) 



|ieA:ieb(0)') 



n 



lieAiiebd)) 



Vi 



(C.69) 



Since A ^ b(0), A has some element in b(0)' Ub(l). Let the largest integer in A be imax- Then, we have Par(E)i^^^^^ - 1, 
and thus, B is odd. This leads to a contradiction. Again, permutations of columns do not affect this discussion. This 
completes the proof of lemma [8] for m-\. 

Appendix C.4. Characteristic vectors 

Let us proceed to the proof for the cases where m > 1. When m - 1, we assigned characteristic values and 1 to 
each column operator according to its symmetry. For m> 1, we will assign a "binary vector" with m components to 
each column operator, which we will call a characteristic vector. We encode "symmetries" of a column operator on 
these characteristic vectors. 

Characteristic vector: We first define two maps /o and /i from !P™^ to for m' > as follows. For a given 
column operator U € !P™'^ which is represented as 



U 



U 



2'"' 



(C.70) 



we define /o(f/),/i(f/)e;P^'„, 

MU)j = UjUj^2«-'-u fi(U)j = Uj 
One may represent fo(U) and fi{U) more explicitly: 



,2"'"'). 



MU) = 



fi(U) = 



U2,..'- 



(C.71) 



(C.72) 



Note that /o and fi decrease the length of the column by half. 

Let us denote a set of all the m component binary vectors as B'"^.. Now, for a column operator U e !P™ ^, we assign 

an m component binary vector b G B^^. through the following rule. 

• If /o(t/) = /, take b,„ = 1 and define f/^'^ = fi{U). 

• IffoiU) + I, take b,„ = and define t/('> = /o(t/). 
and, iterate this procedure: 

• If /o(f/*-'*) = /, b„,-j = 1 and define = /i(t/(^'). 

• If /o(t/<-'') * I, b„,^j = and define t/'^+i^ = foiU^^'')- 

for 1 < j < m - I. We define a characteristic operator V of t/ as follows: 

(m) _ f 4r f iTi\ ^ fnO 



Note that V / when U I. It is worth presenting examples of characteristic vectors here (m = 2): 



(0,0). 



(C.73) 



y 




y 




v 




y 


V 

V 


^ (1,1), 


y 


^ (0,1), 


V 

I 


^ (1,0), 


7 
7 


y. 




/ 




I 




7. 



(C.74) 
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Thus, symmetries of column operators are encoded in characteristic vectors. 
Next, we introduce an order between binary vectors in B'"^.. We define 

m 

8(b) = 2^-2'"' 



(C.75) 



where B is like a binary representation of an integer g(b). For a given pair of m component binary vectors B and b', we 
denote 

B < ^' (C.76) 

if and only if 

8(S) < g{&). (C.77) 

For example, for m = 3, we have the following relations between binary vectors: 

(0,0,0) < (1,0,0) < (0,1,0) < (1,1,0) < (0,0,1) < (1,0,1) < (0,1,1) < (1,1,1). (C.78) 

Below, we shall see that a column operator with a larger characteristic vector is "more symmetric" than a column 
vector with a smaller characteristic vector. 

Property of characteristic vectors: Let us briefly recall the proof of lemma [8] for m - 1 . In the proof, we 
constructed a column operator with b - \ from a column operator with b — O.ln particular, if t/ is a column operator 
with a characteristic value b - Q and a characteristic operator V: 



U 



VV 
V 



where V is some Pauli operator, we have 



UT2(U) 



(C.79) 



(C.80) 



which is a column operator with a characteristic value b = I and a characteristic operator V. Thus, we can create a 
column operator with a larger characteristic value from a column operator with a smaller characteristic value. 

In a way similar to this, one can construct a column operator with b' from a column operator with B as long as 
b' > B. Let us represent a binary column B as follows: 



B = 



B 



1,2'" 



(C.81) 



and denote a set of all the binary columns as B™ Here, we define a parity of B as 

Par{B) = Par{B)i (C.82) 
by viewing a binary column Z? as a binary matrix. For a column operator U € !P™ ^ we define U(B) as follows: 

2'" 



(C.83) 



just like {(B). Note that U{B) is a product of translations of U taken according to a binary column B. Then, the 
following lemma holds. 
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Lemma 9. Consider an arbitrary pair ofm component binary vectors b, b' e B™ ^ such that b < b'. For a given column 
operator U which has a characteristic vector b and a characteristic operator V, there always exists some binary 
column B e B™ , with an even parity Par(B) — such that U (B) has a characteristic vector b' and a characteristic 
operator V. 

In other words, from a column operator U with a characteristic vector b, one can always create a column operator 
U' with larger characteristic vector b' by taking a product of translations of U . On the other hand, it is impossible 
to create a column operator with a smaller characteristic vector from an operator with a larger characteristic vector. 
Therefore, one can create a column operator with higher symmetries (a larger characteristic vector), but cannot 
create a column operator with lower symmetries (a smaller characteristic vector). 

Appendix C.5. Proof of lemma^ 

Now, we prove lemma|9]by explicitly finding a binary column B e B™, for creating U(B) for every pair of b and 
b'. In order to derive such binary matrices, we introduce a certain binary column B{b), called a characteristic column, 
which can be used to change the characteristic vector of a column operator. 

Characteristic column: Given an integer p g Zy", one may have its binary representation by considering the 
inverse of g denoted as : 



P = (Pi 



,Pm) = g \P) 



(C.84) 



where p ^ g(p) ^ YJJ^i Pj2'-'. Now, we define the following sets: 



Js = { 5 e B" 



fl) < bj for all j 



(C.85) 



For example, J(i,o) = {(0, 0), (1,0)) and J(o,o,i) = {(0, 0, 0), (0, 0, 1)). 

Based on Jg, we define the characteristic binary column B{b) e B™ , as follows: 



Bibh.p^i = 



P e if, 
otherwise. 



(C.86) 
(C.87) 



Here, we give some examples: 



J(o.o) = {(0,0)), J(i,o) = {(0,0), (1,0)), 
J(i,i) = {(0,0),(1,0),(0,1),(1,1)) 



'(0,1) 



= {(0,0), (0,1)) 



(C.88) 



and 



B(0,0) = 



B(1,0) = 



B(0,1) = 



B(l,l) = 



(C.89) 



One may see the relation between a characteristic column and a characteristic vector: 



v 




V 




v 




r 


I 
I 


^ (0,0), 


V 

I 


(1,0), 


I 

V 


^ (0,1), 


V 
V 


/ 




I 




/ 




y_ 



(1,1). 



(C.90) 



Thus, if we replace 1 entries in B{b) with V and create a column operator, it has a characteristic vector b. 

In order to discuss changes of characteristic vectors, let us introduce the summation rule between binary vectors. 
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We denote a summation of binary vectors d,b € B"^ asd+b € B™^ and define it as follows: 

g(3) + g(g) = 8(3 + S) (C.91) 

when g(a + b) < 2™ - 1 . Therefore, a + bis just like a summation of two binary "numbers" a and b. 

The characteristic columns defined above can be used to change a characteristic vector of a column operator, as 
summarized in the following lemma. 

Lemma 10. Let b < b'. When U has a characteristic value b and a characteristic operator V, 

U' = UiBiA^)) where S + ^ & (C.92) 

has a characteristic value b' and a characteristic operator V. 

Below, we present a proof of this lemma by finding some property of characteristic columns and a certain rule on 
multiplications of column operators. 

Property of characteristic columns: There is a useful relation between a summation of vectors and characteristic 
columns, as summarized in the following leimna. 

Lemma 11. Let 

2"' 

B{d)*B{$) = YjT2iBid)'^^^^'-jy-^ (mod 2). (C.93) 

7=1 

Then, 

B{d + $) = B{d)*B{$). (C.94) 
The proof involves some exercises on elementary math. 

Proof. We begin by defining a summation of sets of binary vectors. For Bi.Bi, • • • ,B„ c B™^ where a is some 
positive integer, and for b e B™^, consider decompositions: 

^ = + . . . + fo^, fo^ e B, for all i (C.95) 

and denote the number of different decompositions of b as N{b; Bi,B2, ■ ■ ■ ,Bq,). Then, we define the following 
suimnation: 

Bi + B2 + • • • + B„ = { ^ e B^,, : iV(^;Bi,B2, • • • ,B„) = odd }. (C.96) 

With the summation defined above, the claim of the lemma can be written as follows. By setting B = B(d) * B(b), 
we may notice that 

=1 p € Jg + J^, (C.97) 
= otherwise (C.98) 

from a direct calculation. Therefore, we need to show that 

h + h = h.S- (C-99) 

The proof relies on the following sublemma. 
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Sublemma 1. Consider 

aCg = (C.lOO) 

for (2"* > a>p> I). Let the binary representations of a and p be 

S = («!,••• ,a„), $ = (pu--- ,13m) (C.lOl) 

where a — J]™^ a;;2'"' andfi — Xi™ijSi2'"'. Then, aCp is odd if and only if 

Pi < a, for alii. (C.102) 

We suspect that the sublemma above has been proven somewhere else as it seems elementary. Yet, we could not 
find a reference, and thus, we present a proof here. 

Proof. For a given integer p e Z2", let h{p) be the largest integer such that 1^ is an integer. Then, with some 
speculations, one may notice that 

m 

h{p) - ^(2'-i - \)pi (C.103) 

1=1 

where p = {p\,--- ,Pm)- Here, 

a\ 

aCe = , (C.104) 

and 

KaCp) = hia) - h(J3) - hia - P) > 0. (C.105) 

Then, ^Cp is odd if and only if 

h{P) + h{a-P) = h{a). (C.106) 

Let a- p = y. Then, we have 

mm m 

YP-'~' - m + Yj^2'-' - 1)7,- = J](2'-i - l)ai (C.107) 

!=1 1=1 !=1 

and Pi +yi = Oi for all i. This is true if and only if 

Pi < Oi for all i. (C.108) 
This completes the proof of the sublenmia. □ 

Now, let us return to the proof of the lenmia. Below, we prove Ja + = J^^g. Let e{ = (1, 0, • • • ,0). First, we 
show that 

J.1+---+J.1 = h (C.109) 
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where 



S = e1 +--- + e1. (C.llO) 

a 

Here, notice that 

jS G J,^ +---+J,^ (C.lll) 

a 

if and only if ^.C^ is odd since J^^ has two elements: (1, 0, • • • ,0) and (0, • • • , 0). Thus, & J3 from the sublemma, 
and we have 

J.1 +---+J.1 = J^- (C.112) 

a 

Now, let us show that 

J^ + Jg = W (C-113) 
Note that ? e J^^g if and only if a+bCc is odd. Here, we have 



Cc = J]aCi-bCc-i (C.114) 



a+b 

where xCy = when y > x. Let us assume that aC, • bCc-i is odd for / = ai, ■ ■ ■ , ap. Notice that aC, ■ bCc-i is odd 
if and only if both „C, and bCc-i are odd. Then, a+bCc is odd if and only if p is odd. Notice that p is the number of 
decompositions of c such that 

c = C1+C2 (ci G Ja and cl g Jg). (C.115) 
Therefore, p is odd if and only if ? g + Jg. This completes the proof of the lemma. □ 

Multiplication of coluiim operators: Next, let us consider how characteristic vectors change under multiplica- 
tions of column operators. 

Lemma 12. Consider the following column operators: 

U" = UU' i= I (C.116) 

where 

U ^ V, B 

U' V, P . (C.117) 

U" V", b" 



Then, 

• Ifi > &, b" = 9 and V" = V. 

• If 9 =BandV¥= V, b" = Band V" = VV. 

• If 9 =$andV = V, 9' > 9 
Below, we present the proof. 
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Proof. We start with the first claim. When b'^ - 1, b,,, - 1 since b > b' . Then, /o(t/) = fo{U') - /, and fo(U") = /. 
Thus, = 1. When b'„ ^ and b„, = 0, we have fo(U) + I and /o(f/') I. Suppose /()(t/") = /. Then, 
/o(j/) = /o(t/'), and f/(i> = J/'^^'. This means which contradicts with $ > &. Thus, /o(t/") /, and fe,',', = 0. 

When b',„ - and b,n = 1, we have /o(£/) = / and foiU') + /, and fo(U") + /, and b'^ - 0. In summaiy, b'^ - b'„[. One 
can repeat the same discussion and show b' = b" and V" = V. 

The second claim is easy to prove, so we shall skip the proof. Let us move to the third claim. Since 

fhr-- fbJU) ^ fhr-- fbju') = V, (C.i 18) 

we have fh, ■ ■ ■ .//,„,(t/") = /. Let the largest integer / such that fi,. ■ ■ ■ fi,^^^(U") = / be i„ax- Then, ■ ■ ■ fb„,(U") + I 
and b'l -\. If b-,^^^^^ - 1, ■ ■ ■ fbJJJ") - I which leads to a contradiction. Thus, = 0. 
Since ■ ■ ■ fb,„(U") + I, we have 

bi = b'; + l < (C.I 19) 

bi^^^^ = 0, bl^^ = 1. (C.120) 

Thus, If' >b. □ 



Proof of lemma 
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Finally, let us finish the proof of lemma 10 Consider the following column operator V{b) 



which replaces 1 entries in B(b) with V and entries with /: 

V{b)j = y^C'^'J. (C.121) 

Then, one may easily notice that V{b) has a characteristic vector b with a characteristic operator V. Therefore, from 
lemma [TT] we have 

V{$){B{a)) = V{a + ^). (C.122) 
For U with b and V, we decompose U with V{b) as follows: 

U = V(b)U'. (C.123) 



When U + U' has a characteristic vector b' with b' > b from lemma 12 One can repeat the same decomposition 
and obtain: 



where b <b' <b" . Then, from lemma 
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U = V{b)V\b')V'\b")--- (C.124) 
we have 



U{B{a)) = V{d + b)V'{a + b')V"{a + b")--- . (C.125) 



Note that V{c){B{a)) = / if a + c > 2™ - 1. From lemma 12 U{B{a)) is a column operator with a characteristic vector 
d+b and a characteristic operator V. Note that B{a) is an even column when a + (Q, - ■ ■ ,0). This completes the proof 
of lemma [Tol 

Appendix C.6. The existence of an odd matrix for m > 1 

Finally, let us proceed to the proof of lemma [8] to complete the proof of theorem |8] 

Procedure: Consider a logical operator £ defined inside P{x,2"' , 1) with x > 2v. Let us represent characteristic 
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values and characteristic operators for each column of { as follows: 



Uy ^ bu Vy 

^ bl, V,. 



(C.126) 



Without loss of generality, we may assume that fc, < for all i since permutations of columns do not affect parities 
of binary matrices. We define a set of integers such that fc, =bas b(o): 



hib) = [i : bi = b}. 
We denote a group generated by {V,} for i g b(^) as 'Vf 

= {{Vi : i€ hB)}). 
Since x > 2v, there always exists a set of integers A such that 

nv, = /. 



(C.127) 



(C.128) 



(C.129) 



We denote the largest vector in {^,),eA as ba and the largest integer i with Z?, = ba as i = a. Here, we define the 
following binary matrix E{i') 



Eii') = 



such that 



Ei,i = 1 a = n 

Eij = otherwise. 



(C.130) 



(C.131) 
(C.132) 



For example. 



£(1) 



1, 0, 
0, 0, • • • , 

0, 0, • • • , 



£(2) 



0, 1, 0, 
0, 0, 0, 

0, 0, 0, 



(C.133) 



Notice that 

Par{E{i))j = 6ij. 

Now, we define the following operation between a binary column A G ^ and a binary matrix B: 

2" 

B*A = YjTiiB^'-'Y''^ (mod 2). 



(C.134) 



(C.135) 
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Then, consider the following binary matrix: 



E = 2£(0*B(A^) 



(C.136) 



where bj + A/?,- - ba- Note that E{i) * B{Abi) is odd if and only if Abi - (0, ■ ■ ■ ,0). This matrix can generate the 
following column operator: 



J^E(i)*B(Abd 
Y]Ui{B(Abd). 



(C.137) 



ieA 



Note that 



Ui{BiAffi)) 

has a characteristic vector ba and a characteristic operator V, from lemma [ll] Here, we notice that 

Par(E)a = 1 



(C.138) 



(C.139) 



and E is an odd matrix since Aba = (0, ■ ■ ■ , 0). Then, {(E)^- + I since there is no odd identity generating matrix. Now, 
we notice that £{E)x is a column operator with a characteristic vector b'^ > ba from lemma 12 
We summarize the discussion so far as follows. 

• From A such that H/ga - one can form a column operator {(E)^ which has a characteristic vector b'^ > ba 
and a characteristic operator where E is an odd matrix which satisfies 



PariE)a = 1 and Par{E)j = (j > a). 



(C.140) 



Update: Next, we "update" Ua, ba, Va and E(a) to ((E)^, b'^, and E: 



Ua 


^ f(E)x 


bl 




Va 


^ V'a 


E(a) 


E 



(C.141) 



In other words, we replace Ua, ba, Va and E(a) with E(a), {(E)^, b'^, V'a and E, and rename them as t/„, ba, Va and 
E{a). Note that these "updated" E{i) satisfy 

Par{E(i))i = 1 and Par{E(i))j = (; > /)■ (C.142) 

Then, one may repeat the discussion above. Since x > 2v, there is a set A' such that 

Y]Vi = /. (C.143) 

I€A' 

We denote the largest vector in {fe,),eA' as b/j and the largest integer / with bj - bp as / = yS. Then, consider the 
following binary matrix: 

E' = ^£(0*B(A^) (mod 2) (C.144) 

ieA 



62 



where b-, + A/?, = bp. This matrix E' can generate the following column operator: 



Y^E(i)^B(A^d 

ieA' 

Y]Ui{B(Abd). 



(C.145) 



ieA' 



Here, we notice that 

Par{E')p = 1 

since Abp = (0, ■ ■ ■ ,0) and E' is an odd matrix. Then, {{E')jc + I. Note that 

Ui[B{ASd) 



(C.146) 



(C.147) 



has a characteristic vector bp and a characteristic operator V, from lemma[TT| Then, we notice that €{E')^ is a column 



vector with a characteristic vector b'^ > bp. Note that 



Par(E')p = (j>/3). 



(C.148) 



Then, we obtain the following observation. 



From A' such that n/eA" - I for "updated" V,, one can form a column operator ({E')x which has a charac- 
teristic vector b'p > bp and a characteristic operator V^. E' is an odd matrix which satisfies 



Par{E')p = 1 and Par{E')j = {j > p). 



(C.149) 



Here, we again "update" Up, bp, Vp and E(P) to l{E')x, b'^, V'^ and E' . Then, since updated E(i) always satisfy 

Par{E{i))i = 1 and Par{E{i))j = ( j > /) (C.150) 

one can repeat the same discussion again. In each update, characteristic vectors bi increase, and at the end, one ends 
up with the following column operators 



Ui 1, Vi, £(1) 

t/2 ^ 1, Vi, E(2) 

U, f, y,, E(x) 

where 1 s (1, ■ ■ ■ ,1) and 

Par{E{i))i = 1 and Par(E(i))j = (; > 
Then, there exists a set A such that 

ISA 



(C.151) 



(C.152) 



(C.153) 
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and, the following matrix is an identity generating matrix 



E = Y^Eii) (mod 2). (C.154) 

Let the largest integer in A be i,nax- Then, E is odd since Par{E)i^^^^^ - 1. However, this contradicts with our original 
assumption that there is no odd identity generating matrix. This completes the proof of lemma |8] lemma [6] and 
theorem [8] 

Appendix D. Derivation of logical operators 

Having showed that two-dimensional logical operators can be decomposed as a product of two-dimensional and 
one-dimensional centralizer operators, let us proceed to the proof of theorem |5] The proof owes a lot to arguments 
presented in IJ7J. For simplicity of presentation and in order to avoid making the paper unnecessarily long, we shall 
skip some parts of the derivation. However, we believe that interested readers can easily construct rigorous proofs. 

Preliminaries: We begin by providing some corollaries and lemma which are useful in the derivations of logical 
operators. Let us first generalize theorem|8]for any nj. 

Corollary 2. Consider a three-dimensional STS model with the system size n\ =2-2^"'^'!, arbitrary n2 and n^, > 1. For 
a given logical operator £ supported inside P(ni,n2, 1), one can decompose £ as a product of the following centralizer 
operators 

( ~ tatb, {a,fh e C/>(„, ,„,,!) (D.l) 

where 

= 4, where /? < 2^"^" (D.2) 

and £a is defined inside P(2v, n2,l). 

The proof relies on the fact that one can make a logical operator "quasi-periodic". 

Proof. Let us represent n2 as «2 - 2"' ■ n'j where n'^ is some odd integer Then, for a logical operator I defined inside 
P(ni,n2, 1), one can see that the following logical operator 

'A 

{' = [^r*-''""^"(/') ~ ( (D.3) 

7=1 

is equivalent to ( since {' is a product of an odd number of translations of £. Notice that {' is periodic in the 2 direction: 

rf (f) = e (D.4) 

with the periodicity 2™. Because of this periodicity, one can form identity generating matrices in a way similar to the 
cases when n2 - 2". Therefore, one can see that theorem [8] holds for any n2- Q 

We have seen that a two-dimensional logical operator can be decomposed as a product of a one-dimensional cen- 
tralizer operator and a two-dimensional centralizer operator, as summarized in theoremjS] One can further decompose 
a one-dimensional logical operator as a product of a one-dimensional centralizer operator and a zero-dimensional 
centralizer operator, as summarized in the following lemma. 

Lemma 13. Consider a three-dimensional STS model with the system size ni = 2 ■ 2'^"'-''\, n2 = 2 ■ 2^^'^'\ and 
«3 > 1 where m is an arbitrary positive integer. For a given logical operator I supported inside P(2v, n2, 1), one can 
decompose { as a product of the following centralizer operators 

i ~ (a(b^ fa^fh 6 Cp(2v,n2,l) (D-5) 
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where 



where p < 2<2>')' 



(D.6) 



and €a is defined inside P(2v, (2v) , 1). 



We show the claim of lemma [T3lg raphically in Fig. |D.34| One can prove the lemma through discussion similar to 
the one used in the proof of lemrnajs] So, we shall skip the proof 



(a) 



1^ ly 



Pi2v,n2,l] 



f-1 — rn ^ 



















la 







{2v) 



periodic in 2 direction 



2v 



Figure D.34: The claim of lemma [T3] 

Finally, let us extend the claim of theorem [4] slightly. 
Corollary 3. Consider a two-dimensional STS model with even n\ and n^. 
• Let ( be a logical operator which is periodic in the 1 and 2 directions: 



(D.7) 



Then, { is a two-dimensional logical operator, and there exists a zero-dimensional logical operator r which is 
defined inside 2v) and anti-commutes with {. 



Let { be a logical operator which is defined inside P(\,n2) and periodic in the 2 direction: 

T2({) = £. 



(D.8) 



Then, t is a one-dimensional logical operator, and there exists another one-dimensional logical operator r 
which is defined inside P{n\, 1), anti-commutes with {: {(, r} — and periodic in the 1 direction: 



Ti(r) 



(D.9) 



In other words, if we find a logical operator which is periodic in both 1 and 2 directions, we readily know that 
it is a two-dimensional logical operator. Similarly, if we find a logical operator which is defined inside P{l,n2) and 
periodic in the 2 direction, we readily know that it is a one-dimensional logical operator Note that the corollary holds 
only for the cases where both ni and are even. 
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Proof. Since £ is periodic and system sizes are even, { commutes with all the one-dimensional logical operators and aU 
the two-dimensional logical operators. Therefore, ( can anti-commute only with zero-dimensional logical operators. 
This means that { is a two-dimensional logical operator The second claim can be proven in a similar way. □ 

Now, we derive all the logical operators in a three-dimensional STS models, as described in theorem|5] Consider 



the system size analyzed in lemma 13 (ni, n2,n3) = (2 • 22"2>'!,2 • 2(2''> !,«3) where «3 > 1 is some fixed integer. 



Here, we view the entire system as a two-dimensional system by considering f (1, n2, 1) as a single composite particle 



(Fig. D.35 I. (So, the entire system is viewed as a two-dimensional lattice of one-dimensional tubes). 



For a two-dimensional STS model, we already know the geometric shapes of all the logical operators, as sum- 
marized in theorem |4] When viewed as a two-dimensional system, "zero-dimensional" logical operators are defined 
inside P{2vn2,n2, 1). From theorem |8] we notice that these "zero-dimensional" logical operators can be actually 
defined inside Pi2v, «2j !)■ Then, we have the following logical operators (Fig. D.35 1. 



(a) 



(b) 



Figure D.35: Viewed as a two-dimensional system. 



An anti-commuting pair of logical operators in Fig. |D.35{ a) where i is defined inside P(2v,n2, 1) and r is 
periodic: 

Tdr) = n{r) = r. (D.IO) 
An anti-commuting pair of logical operators in Fig. |D.35| ^b) where £ is defined inside f (1, «2, ns) and periodic: 



and r is defined inside P{n 1,112, 1) and periodic: 



(D.ll) 



(D.12) 



Below, we analyze anti-commuting pairs of logical operators in Fig. D.35 a) and Fig. D.35 b), and derive logical 
operators. 

Pairs in (a): Below, we analyze properties of logical operators described above. We start with anti-commuting 
pairs described in Fig. D.35 a). We stop viewing the system as a two-dimensional system for the moment. Logical 
operators defined inside P(2v, 112, 1) consists of periodic one-dimensional logical operators and zero-dimensional log- 
ical operators defined inside P{2v, (2v)^, 1) due to lemma 13 Let us first analyze a zero-dimensional logical operator 
( defined inside P(2v, {2vf, 1). 
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From theorem|4] ( is also a logical operators for arbitrary n\ and «3. Consider the case when n-i - \ (Fig. D.36| l. 
Then, by viewing the system as a two-dimensional system which extends only in the 1 and 2 directions, one notices 
that there exists a two-dimensional logical operator r which is periodic in both 1 and 2 directions: 

Tx(r) = T2(r) = r (D.13) 

and anti-commutes with £: {t,r] - 0. Next, let us consider the case when > 1. Then, one may extend the 



construction of r as follows (Fig. D.36 1: 



(D.14) 



x=l 



In other words, we put r in a periodic way in the 3 direction to form r' . We shall call such an extension the periodic 
extension. The three-dimensional logical operator r' obtained after the periodic extension of r is periodic in all the 
directions: 

Tdr') = T2(r') = T^ir') = / (D.15) 
and anti-commutes with { : {r', £] - 0. Then, one may notice that r' and { form a pair of anti-commuting logical 



operators for any system size n. From this discussion, we obtain the following observation (Fig. D.36 1 



For zero-dimensional logical operators defined inside P{2v, (2v)^, 1), there always exists a three-dimensional 
logical operator r which is periodic: 



Tiir) = T2{r) = T^{r) = r 
and anti-commutes with i: {(, r] -Q. ( and r are logical operators for any system size n. 



(D.16) 
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Figure D.36: Constructions of zero-dimensional and tliree-dimensional logical operators. 

Next, let us consider a one-dimensional logical operator I defined inside P{2v,n2, 1) which is periodic in the 2 
direction: 



(D.17) 



67 



for j6 < 2*^'''\ Recall that r is periodic: 

T,{r) = W = r (D.18) 

and anti-commutes with £ : {r, {} - 0. Then, one can decompose r as a product of two centralizer operators from 
lemma |5] 

r ~ r„rh (D.19) 

where 

Tiira) = TsCrJ = r«, TiCr^) = TsCr^) = (D.20) 
and is defined inside P(ni,2v,nj,), and 

(r,) = r, (D.21) 

for jS' < 2^^'. Then, we notice that 

[rbJ] = (D.22) 

since T2({) - t and (rt,) - rt, and nxIPP' is an even integer for n2 = 2 ■ 2^^'-'' !. Thus, we have 

{raJ} = 0. (D.23) 

Since is periodic in the 1 and 3 directions, one can periodically extend its construction for arbitrary ni and 
M3. Now, let us consider the system size such that n\ and n-i are odd. Here, note that has some equivalent logical 
operator r'^ defined inside P(n 1 , 1 , ns) 1 17 1 . Then, the following operator 

C = \\T[TiiO ~ (D.24) 

•J 



is equivalent to (Fig. D.37 1. Note that r'^ is defined inside P(ni, 1, ns) and periodic in the 1 and 3 directions: 

TdO - TiiO = r'^. (D.25) 

Finally, we show that £ can be periodic in the 2 direction. Consider the case when «i and «3 are even, and n2 - 1. 
Then, r'J is also a logical operators. Now, there always exists some logical operator £' defined inside f (1, 2v', 1) which 
anti-commutes with r'^ from corollary |3] One can periodically extend its construction to arbitrary n2- We denote it i" . 
Then, {" and r" are logical operators for any system size. From this discussion, we obtain the following observation 
(Fig. [DJ7] i. 

• For one-dimensional logical operators defined inside P(2v, 112,!), there exists a two-dimensional logical opera- 
tor r defined inside P(ni , 1 , n^) which is periodic: 

Tdr) = T,(r) = r (D.26) 

and anti-commutes with f: {{, r} = Q. { can be also periodic: 

Tiit) = £, (D.27) 

and, { and r are logical operators regardless of the system size n. 



Pairs in (b): Let us proceed to the analysis on pairs of logical operators in Fig. D.35 b). We consider the following 
anti -commuting logical operators £ and r. 
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Figure D.37: Constructions of one-dimensional and two-dimensional logical operators. 



• lv& defined inside P(l,n2, n^) and periodic: T^ii) - (. 

• r is defined inside P{n\, ni, 1) and periodic: T\{r) - r. 

Since ( is periodic in the 3 direction, one can decompose ( as follows: 

i ~ 44, T^ifa) = 4 and 73(4) = 4 
where 4 is defined inside f (1, 2v, ^3), and 4 is defined inside P(l,«2, 13) and periodic: 

rf(4) = 4 



(D.28) 



(D.29) 



where fi < 2^''. Thus, logical operators defined inside P(l,«2,n3) consist of two-dimensional logical operators and 
one-dimensional logical operators defined inside f (1, 2v, n-i,). 

Let us analyze a one-dimensional logical operator ( defined inside P(l,2v, «3) first. We decompose r defined 
inside P(ni, «2, 1) as follows: 



r ~ ran, 

where is defined inside P(«i,2v, 1), and rb is defined inside P(ni,«2, 1) and periodic: 

T2(rb) = rh 

where /3' < 2^''. Then, we notice that 

[{,ra] = 

since there exists a translation of which does not overlap with {. Thus, we have 

14^^,1 = 0. 



(D.30) 



(D.31) 



(D.32) 



(D.33) 



Let us consider the case where ni = 1, and «3 is even. Note that C and r/, are both logical operators. Then, from 
corollary [3] there exists a one-dimensional logical operator r' which is defined inside P{l,n2, 1), anti-commutes with 
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I and is periodic in the 2 direction: 

r2(/) = /. (D.34) 

Then, we periodically extend r' in the 1 direction and define r" . Then, we notice that r" is defined inside P(ni,n2, 1) 
and periodic in the 1 and 2 directions. From this discussion, we obtain the following observation. 

• For a one-dimensional logical operator { defined inside P(l,2v, n^), there exists a two-dimensional logical 
operator r defined inside P(n\,n2, 1) which is periodic: 

Tiir) = T2(r) = r (D.35) 

and anti-commutes with (: {(, r] -Q. ( can be also periodic: 

Tii{) = (, (D.36) 

and, I and r are logical operators regardless of the system size n. 

Finally, let us analyze a two-dimensional logical operator t defined inside P(\,n2, n^) with 

rf(^) = (D.37) 

Since r is periodic in the 1 direction, one can decompose it as follows: 

~ r^rh, Ti(ra) = r„ and Ti(rh) = r^, (D.38) 

where ra is defined inside P{ni , 2v, 1) and rt, are defined inside P(ni , n2, 1) and periodic: 

(r,) = (D.39) 

where fi' < 2^''. Then, one may notice that 

{£,rb} = 0. (D.40) 
Then, the rest is immediate, and we obtain the following observation. 

• For a two-dimensional logical operator i defined inside P{l,n2,nj,), there exists a one-dimensional logical 
operator r defined inside P(ni, 2v, 1) which is periodic: 

Ti(r) = r (D.41) 

and anti-commutes with {: {{, r} = 0. f can be also periodic: 

T2(() = T,(£) = £, (D.42) 

and, { and r are logical operators regardless of the system size n. 



Let us recall the discussion so far. We started our analysis with the system size considered in lemma 13 Then, for 
each pair of anti-commuting logical operators, we found logical operators whose geometric shapes are the same as 
the ones described as in theorem [5] These logical operators are also logical operators for other system sizes because 
of their periodic structures and scale symmetries of the system. 

It remains to sort these logical operators in a canonical form by analyzing their commutation relations. However, 
we shall skip this process since it is straightforward. See L17J for a similar discussion. This completes the proof of 
theorem |5] 
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